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Synchrotron signature of a relativistic blast wave with decaying 
microturbulence 
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ABSTRACT 

Microphysics of weakly magnetized relativistic collisionless shock waves, corroborated by 
recent high performance numerical simulations, indicate the presence of a microturbulent 
layer of large magnetic field strength behind the shock front, which must decay beyond some 
hundreds of skin depths. The present paper discusses the dynamics of such microturbulence, 
borrowing from these same numerical simulations, and calculates the synchrotron signature 
of a powerlaw of shock accelerated particles. The decaying microturbulent layer is found to 
leave distinct signatures in the spectro-temporal evolution of the spectrum F v oc t~ a v~P of 
a decelerating blast wave, which are potentially visible in early multi-wavelength follow-up 
observations of gamma-ray bursts. This paper also discusses the influence of the evolving 
microturbulence on the acceleration process, with particular emphasis on the maximal energy 
of synchrotron afterglow photons, which falls in the GeV range for standard gamma-ray burst 
parameters. Finally, this paper argues that the evolving microturbulence plays a key role in 
shaping the spectra of recently observed gamma-ray bursts with extended GeV emission, such 
as GRB090510. 
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1 INTRODUCTION 

The acceleration of particles at a decelerating relativistic collision- 
less shock front constitutes a key building block of the afterglow 
model of gamma-ray bursts (GRB, Meszaros & Rees 1997). The 
standard phenomenology models the accelerated electron popula- 
tion as powerlaw d/V e /d7 e oc 7,T P , which radiates powerlaw pho- 



ton spectra of the form F v oc t 



with a temporal decay in- 



dex a and a frequency index /3 that are direct functions of p, see 
e.g. Piran (2005) for a review, or e.g. Sari et al. (1998), Panaitescu 
& Kumar (2000) for detailed formulae. From both microscopic and 
observational points of view, the situation however appears more 
complex, in spite of several remarkable results of the past decade. 

On the microscopic level, for instance, one understands the 
formation of a relativistic collisionless shock front in a weakly 
magnetized medium - such as the interstellar medium (ISM) - 
through the self-generation of intense small scale electromagnetic 
fields that act as the mediating agents for the transition from the 
far upstream unshocked state to the far downstream shocked state. 
The accelerated particles, as forerunners of the shock front, play a 
central role in triggering the microinstabilities that build the self- 
generated field. In turn, this self-generated microturbulence con- 
trols the scattering of these particles and it therefore directs the ac- 
celeration process, which becomes intricately non-linear. This gen- 
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eral scheme has been validated so far in high performance particle- 
in-cell (PIC) simulations (e.g. Spitkovsky 2008, Keshet et al. 2009, 
Sironi & Spitkovsky 2009, 2011, Martins et al. 2009) and under- 
stood on the basis of analytical arguments at the linear level (e.g. 
Medvedev & Loeb 1999, Gruzinov & Waxman 1999, Lyubarsky 
& Eichler 2006, Lemoine & Pelletier 2010, 2011a, Rabinak et 
al. 2011). The situation becomes more complex when one tries 
to bridge the gap between the limited simulation timescales and 
the much longer timescales probed by the observations. On such 
timescales, one would indeed expect that the microturbulence has 
died away (e.g. Gruzinov & Waxman 1999), yet GRB afterglow 
modelling has generally pointed to a field strength close to a per- 
cent of equipartition permeating the blast, on day timescales. The 
origin of this field and its relation with the microturbulence behind 
the shock front has remained a nagging issue for many years. 

On the observational level, the recent era of rapid follow- 
up observations in the X-ray and GeV domain has brought its 
wealth of surprises. The Swift satellite has revealed X-ray after- 
glow light curves that differ appreciably in the 10 2 — 10 4 s do- 
main from the canonical afterglow model (Nousek et al. 2006, 
O'Brien et al. 2006). Of more direct interest to the present work, the 
Fermi-LAF telescope has reported the discovery of long-lived (~ 
100—1000 s) GeV emission in a fraction of observed bursts. In one 
case (GRB090510), this emission has been measured almost con- 
temporarily to emission in the X-ray and optical domains as early 
as 100 s (Ackermann et al. 2010, de Pasquale et al. 2010). This 



2 M. Lemoine 



long-lived high energy emission has been shown to fit nicely the 
predictions of a model in which the electrons cool slowly through 
synchrotron radiation in a background shock compressed magnetic 
field, without any need for microturbulence (Barniol-Duran & Ku- 
mar 2009, 2010, 2011a; see also Gao et al. 2009; de Pasquale et 
al. 2010; Corsi et al. 2010; Ghirlanda et al., 2010; He et al 2011; 
and see Ghisellini et al. 2010, Razzaque 2010 and Panaitescu 201 1 
for alternative points of view). Given the past history in GRB after- 
glow modelling, such a low magnetization of the blast may come 
as a surprise, but it may also point out that after all, the microturbu- 
lence does decay away as theoretically expected, and that the high 
level of turbulence seen on day timescales has been seeded through 
some other instabilitjQ. 

Depending on how fast and how far from the shock this micro- 
turbulent layer decays, it is likely to influence the particle energy 
gains from Fermi acceleration, and losses through synchrotron ra- 
diation. The microturbulent layer must actually ensure the scatter- 
ing of accelerated particles, because in the absence of microturbu- 
lence, these particles would be advected away with the transverse 
magnetic field lines to which they are tied and acceleration would 
not take place (e.g. Begelman & Kirk 1990, Lemoine et al. 2006, 
Niemiec et al. 2006, Pelletier et al. 2009fl Now, scattering in small 
scale turbulence is so slow that producing GeV photons at an ex- 
ternal blast wave of Lorentz factor of a few hundreds represents 
a challenge, see e.g. Kirk & Reville (2010), Lemoine & Pelletier 
(201 lc), and see also Piran & Nakar (2010), Sagi & Nakar (2012). 
It would be about impossible if the particles were to scatter in the 
microturbulent layer then radiate in the much weaker shock com- 
pressed background field. Therefore, the interpretation of Barniol- 
Duran & Kumar (2009, 2010, 201 la) actually suggests that the mi- 
croturbulence also plays a role in the radiation of GeV photons, at 
the very least, if not in shaping the synchrotron spectra over the 
broad spectral range. In other words, the observation of extended 
GeV emission and its early follow-up in other wavebands may be 
opening a rare window on the dynamics of the microturbulence in 
weakly magnetized relativistic collisionless shocks. 

In this context, the present paper proposes to discuss the 
synchrotron spectra and more generally the afterglow spectrum 
F u oc t a v~P of a decelerating relativistic blast wave, account- 
ing for the time evolution of the microturbulence behind the shock 
front. While the initial motivation of this work was to provide a 
concrete basis for the scenario of Barniol-Duran & Kumar (2009, 
2010, 2011a), in which particles scatter in a time decaying micro- 
turbulence but radiate in a region devoid of microturbulence, it has 
become apparent that the possibility of decaying microturbulence 
opens a quite rich diversity of phenomena, which deserve to be dis- 
cussed in detail. This problem has been tackled by Rossi & Rees 
(2003), who considered the simplified case of a homogeneous mi- 
croturbulent layer that dies instantaneously beyond some distance, 
and by Derishev (2007), who showed that a particle radiating in a 
time evolving magnetic field can lead to spectra quite different from 
the standard one-particle synchrotron spectra. Borrowing from the 
latest PIC simulations, the present paper establishes a model of the 
microturbulence strength and coherence length evolving as power- 
laws in time beyond some distance, until the decay saturates down 

1 De Pasquale et al. 2010 and Corsi et al. 2010 have shown that the 
afterglow emission could be modelled with a more traditional estimate 
eg ~ 10~ 2 — 10~ 3 , but this comes at the price of an extraordinarily 
low external density n < 10 — 6 cm -3 . This interpretation is not considered 
here, see also Sec. |3.3| for further discussion. 
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to the background shock compressed magnetic field; the present 
study then calculates the synchrotron spectra in this structure for 
various typical configurations (slow cooling, fast cooling, with and 
without inverse Compton losses) and it discusses the problem of 
particle scattering, acceleration timescale and maximum photon en- 
ergy in this setting. As such, it generalizes and encompasses these 
former studies, in the spirit of providing new tools with which one 
can analyse existing and forthcoming data. A brief comparison to 
present early afterglow observations is provided. 

The detailed spectra and the spectro-temporal indices a, j3 of 
F v oc t~ a v~^ are provided in Appendix lAl while Section [2] de- 
tails the model for the evolution of microturbulence and provides 
the general characteristics of the afterglow light curves and spectral 
energy distributions in various configurations. Section |3]discusses 
the scattering process and the maximal acceleration energy, and it 
confronts the above models to existing data. The results are summa- 
rized in Sec. [4] Throughout, this paper adopts the standard notation 
Q x = Q/10 x with Q a generic quantity in cgs units. Fiducial val- 
ues used for numerical applications correspond to those derived in 
Barniol-Duran & Kumar (2009, 2010, 201 la), and He et al. (201 1), 
e.g. an external density n ~ 10~ 3 cm~ 3 , a blast Lorentz factor 
7b ~ 300 at 100 s and it is assumed that the blast has entered 
the deceleration regime. One must distinguish the time t bs in the 
observer frame, also written t ^ B = 100 ti s in numerical applica- 
tions, from the time experienced by a particle since shock entry; 
this difference is manifest everywhere. For convenience, the nota- 
tion z+,o.3 = (1 + • z )/2 is introduced, z denoting the redshift of 
the GRB. 



2 SYNCHROTRON SPECTRA WITH TIME DECAYING 
MICROTURBULENCE 

The self-generation of microturbulence in the precursor of a rel- 
ativistic collisionless shock front propagating in a very weakly 
magnetized medium appears both guaranteed and necessary to the 
maintenance of the shock. It is necessary because in the absence 
of a background magnetic field, self -magnetization is required to 
build up a magnetic barrier that initiates the shock transition. It 
is guaranteed because the development of microinstabilities fol- 
lows naturally from the penetration of the unshocked plasma by 
the anisotropic beam of supra-thermal particles moving ahead of 
the shock (e.g. Medvedev & Loeb 1999). 

Studies of non-relativistic collisionless magnetospheric 
shocks have shown that dissipation ahead of the shock is initi- 
ated by the reflection of ambient particles (i.e. from the unshocked 
plasma) on the shock front in the compressed magnetic field (e.g. 
Leroy et al. 1982). Particle-in-cell simulations indicate that a sim- 
ilar phenomenon takes place at a relativistic unmagnetized shock, 
although the particles now reflect on the small scale electromag- 
netic fields self-generated by the microinstabilities (e.g. Spitkovsky 
2008). The reflected and accelerated particle populations merge to- 
gether and trigger microinstabilities such as the Weibel (fi lamen- 
tation) instability or oblique electrostatic instabilities, provided the 
precursor extends far enough for these modes to grow on the pre- 
cursor crossing timescale (Lemoine & Pelletier 2010, 2011a). In a 
very weakly magnetized shock wave, with magnetization typical of 
the ISM and blast Lorentz factor 7b < 10 3 , this appears guaran- 
teed. 

As seen from the shock frame (in which the shock front lies 
at rest) the incoming kinetic energy is carried by the protons, 
the electrons carrying only a fraction m e /m p of the incoming 
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flow kinetic energy. Energy transfer between the two species in 
the microturbulence leads to heating of the electron population, 
close to equipartition by the time it reaches the shock front, as 
observed in current PIC simulations (Sironi & Spitkovsky 2011), 
see also the discussion in Lemoine & Pelletier (2011a). Equipar- 
tition means that the incoming electrons carry Lorentz factor 
7e ~ 7bm p /m e , hence their skin depth scale (downstream frame) 
c/ [47T7bn u / (7 e m e )] 1//2 ~ c/oj p j. The natural length scale of the 
electromagnetic structures produced by these microinstabilities is 
therefore the ion skin depth scale c/oj p i of the upstream plasma. 



2.1 Input from Particle-in-cell simulations 

Particle-in-cell simulations not only validate the above general 
scheme, they also provide interesting constraints on the shape and 
evolution of microturbulence ahead and behind the shock front. 
Two most recent and most detailed studies are of direct interest 
to the present work. 

Chang et al. (2008) have performed long simulations of the 
evolution of micro-turbulence behind a relativistic shock front. The 
simulations have been computed for an unmagnetized pair plasma 
with relative Lorentz factor 7b = 15 between upstream and down- 
stream. In weakly magnetized relativistic shock waves, the preheat- 
ing of electrons in an electron-ion shock of similar configuration 
appears so efficient that for all practical matters, the downstream 
plasma behaves as a relativistic pair plasma; hence the results of 
Chang et al. (2008) can be transposed to an electron-ion shock. 
These simulations show that the microturbulence remains mostly 
static in the downstream rest frame, and that it is composed of an 
intermittent magnetic field structure that can be roughly described 
as a collection of magnetic loops and islands on typical length 
scales ~ 10 — 30 c/cj p j. One clear observation made in this work is 
that the small scale structures dissipate first, leaving the large scale 
clumps unaffected over the timescale of the simulation. Chang et al. 
(2008) interpret this gradual erosion as collisionless damping, with 
a damping frequency Qui oc A -3 (A denoting the spatial scale). If 
the magnetic turbulence is described in Fourier space as a power 
law spectrum with most of magnetic power on small length scales, 
this implies a decay of the magnetic field strength accompanied by 
an evolution of the coherence scale Xsb oc t 1 ^ 3 . This is made ex- 
plicit further below. 

The longest PIC simulation so far for a relativistic shock is 
that of Keshet et al. (2009), which extends to about lO 4 ^"; 1 ~ 

— 1/2 

240 w_ 3 s (comoving time). For values envisaged by Barniol- 
Duran & Kumar (2009, 2010, 2011a) and He et al. (2011) to de- 
scribe the GeV extended emission, i.e. an ejecta of energy E ~ 
10 53 erg, launched into a medium of density n ~ 10~ 3 cm~ 3 with 
initial Lorentz factor 7 c j ~ 10 3 , the above timescale represents 
close to 1 % of a dynamical timescale at an observer time of a hun- 
dred seconds. Keshet et al. (2009) provide a detailed study of the 
magnetic field power spectrum of the turbulence and its evolution. 
They confirm most of the findings of Chang et al. (2008); in partic- 
ular, they show that the magnetic field does decay behind the shock 
wave, but on rather long length scales compared to a skin depth 
c/ojpi. More importantly, they find that the presence of shock ac- 
celerated particles influences the decay timescale of the magnetic 
field with a general trend being that higher energy particles ensure a 
longer lifetime for the downstream microturbulence. Given that the 
simulations of Keshet et al. (2009) extend for a time that is much 
smaller than the dynamical time, it has not had time to produce very 
high energy particles and to probe their impact. Such high energy 



particles would tend to populate the magnetic perturbation spec- 
trum with longer wavelength modes, which would then decay on 
longer timescales when downstream. In any event, this should not 
call into question what has been said above, since low energy par- 
ticles carry most of the energy of a shock accelerated population 
with index p > 2. Nevertheless, to probe how far the perturbation 
spectrum may be populated, one can conduct the following exer- 
cise. The maximal size of the precursor is given in a reasonable 
approximation by r L ,max,i,o/7sh> where r L ,max,i,o represents the 
gyration radius of the highest energy ions in the background up- 
stream magnetic field, and 7 s h = V^Jb the shock Lorentz factor 
as measured upstream. This result is discussed in detail in Plot- 
nikov et al. (2012) but it can be understood as follows. The high- 
est energy ions are those that travel the furthest away from the 
shock, since electrons are generally accelerated to a smaller en- 
ergy due to synchrotron losses; furthermore, the particles gyrate by 
an angle l/7 s h over a timescale t rcs , u ~ c _1 r- L , max, j,o/7sh be- 
fore being caught up by the shock front (Achterberg et al. 2001); 
finally, the typical distance between the shock front and the par- 
ticle is ct rcs , u (1 — p sh ) ~ c£ rcSi u/(27s h )- Assuming that the 
ions are accelerated on a timescale t rcSiU (see also Sec. I3.lt and 
comparing the time available for acceleration with the age of the 
shock wave r/c, one finds the precursor size ~ r/j^ h ~ 2 x 
10 4 (c/wpOnVs 2 i2 2_i_o 3- This indicates that the perturbation may 
well extend on several decades. In the absence of ions, the precursor 
size would be set by the highest energy electrons; balancing accel- 
eration at a Bohm rate (as experienced downstream) to synchrotron 
losses, one would find a precursor size about 20 times smaller for 
the adopted fiducial parameters, nevertheless much larger than the 
typical size of the fluctuations. 

Finally, Keshet et al. (2009) suggest a damping frequency 
Qlu oc A~ 2 , which implies, when combined with the above result 
that the decay of the microturbulence might extend over quite long 
spatial scales. 

Following Chang et al. (2008), the magnetic field power spec- 
trum is described as a time decaying powerlaw form in the down- 
stream (comoving) frame, with 

(SB(tf) =a B 5Bl £ maX g(Ay B exp (_£) , (1) 

with t denoting the time (downstream frame) since shock entry of 
the corresponding plasma element, A M (resp. A max ) the minimum 
(resp. maximum) wavelength scale of the microturbulence at t = 0, 
as < —1 so that the turbulent power lies at the smallest scales, 
as = 1 1 + qlb I for normalization purposes, SB^ denotes the rms 
field strength at t = and t\ — jSscjj -1 the damping time, which 
depends on A: 

t x = Up; 1 (oj pi A/c) aA . (2) 

Assuming A max 3> A M for the moment, Eq. (QJ can be inte- 
grated in terms of an incomplete Gamma function, 

(5B(tf) = ^ QA -V(t) (1+ ° B)/ai 

f / l + a B \ [ 1 + a.B ,.-11 
x < T r ; fi(t) \ , 

(3) 

with 

/i(t) = u pi t ^ • (4) 

For convenience, one may approximate Eq. ((3} with respectively 
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the small and large argument limits to describe the evolution of the 
magnetic field strength as 



(SB(t) 2 ) 



SBl if n{t) < 1 , 

5BlT (1 + \a t \) ti(t) at if M (t)»l, 



with the following definition 



a, 



1 + CtB 
a\ 



< . 



(5) 



(6) 



In the following, the numerical factor V (1 + \at\), of order unity, 
will be dropped henceforth. 

Equation [6] shows that the temporal decay index of the mag- 
netic field behind the shock is inherently linked to how power is 
distributed on scales larger than the minimum scale A M - as charac- 
terized by olb - and to how fast small scale features are dissipated 
- as characterized by a\. The interpretation for this is clear: as 
small scales are erased, magnetic power is removed, but the rate at 
which the total strength erodes depends on how much strength is 
left at longer wavelengths. In fine, the uncertainty on a t is related 
to the sourcing of large wavelength fluctuations, which are likely 
related to the dynamics of high energy particles in the upstream. 
The above assumption A max 2> A M has been discussed above. Its 
robustness depends crucially on the influence of maximal energy 
particles upstream of the shock front. In the extreme opposite case 
A m ax ~ A M , one should observe a roughly constant magnetic field 
while fj,(t) < 1, followed by fast decay once fi(t) > 1. This situa- 
tion may be accounted for by Eq.[5]wifh a more pronounced value 
of at- 

The following therefore considers a range of possibilities for 
at, even though the PIC simulations of Chang et al, (2008) and 
Keshet et al. (2009) both suggest — 1 < a t < 0. More specifi- 
cally, Chang et al. (2008) suggest that the magnetic field Fourier 
spectrum (for SB, not SB 2 ) in wavenumber has slope ~ — > 1/2, 
which corresponds to 1 + as ~ —2 — > — 1, and a\ = 3 leading to 
at ~ -1/3 -» -2/3. Keshet et al. (2009) show that right behind 
the shock front, the magnetic field decays exponentially on short a 
distance scale to level off at a strength corresponding to eb ~ 10~ 2 
for some hundreds of skin depttjf]. A closer inspection of their Fig. 3 
however reveals that the initial exponential decay leaves way to a 
powerlaw decay at late simulation times (thus meaning far down- 
stream) and by eye, one estimates at ~ —0.5. These simulations 
thus indicate a value of at between —1 and 0; however, given the 
present limitations of the PIC simulations, and the above possible 
caveat related to the extension of the magnetic perturbation spec- 
trum, one cannot exclude yet that at < —1. In this respect, one 
must point out that recent simulations of the development and the 
dynamics of relativistic Weibel turbulence indeed suggest a value 
at — —2 (Medvedev et al. 2010). Although these simulations do 
not simulate the shock itself, but a Weibel turbulence through the 
interpenetration of two relativistic beams, these are 3D while the 
shock simulations of Chang et al. (2008) and Keshet et al. (2009) 
are 2D. 

In order to account for these different possibilities in the fol- 
lowing, Eq.[5]is kept in its present form, but the decay exponent at 
is assumed to take possibly mild or more pronounced values. De- 
pending on whether at < — 1 or — 1 < a t < 0, it will be seen that 
radically different radiative signatures are to be expected. 

3 This result motivates the present choice of eg = 10~ 2 as a fiducial 
value, even though the magnetic energy density reaches ~ 15 % of the 
incoming energy at the shock transition itself (see Chang et al. 2008, Keshet 
et al. 2009). 



In summary, theoretical arguments combined with recent high 
performance PIC simulations suggest the following characteri- 
zation for the evolution of the microturbulence behind a rela- 
tivistic (weakly magnetized) shock front. Immediately behind the 
shock, the magnetic field carries strength SB^ corresponding to an 
equipartition parameter es = SB 2 ./ (32-Ky 2 nm p c 2 J with fiducial 
value es ~ 10~ 2 , while A M ~ 10 — 30c/tJ p i represents the fiducial 
value for the coherence scale at that same location. The magnetic 
field strength decays as t at after a time 



i M+ = cjpi 1 (^ pi A M /c) c 



(7) 



defined through p(t M _|_) = 1, of the order of hundreds to thou- 
sands of inverse plasma times, until it eventually settles at the shock 
compressed value B^ = 47b-B u . In the following, this timescale is 
rewritten in units of inverse plasma times as 



(8) 



meaning also that the undecayed part of the microturbulence ex- 
tends for A M skin depths. Note that A M 3> 1 according to the 
above simulations. Finally, the coherence length of the microturbu- 
lence evolves as t 1 ' ax , with fiducial value ~ 2 — 3. 



2.2 Radiation in time decaying microturbulence 

As a particle gets Fermi accelerated, it interacts with the turbulent 
layer within a scattering length scale i sca tt of the shock front. This 
scattering length scale controls the residence time hence the accel- 
eration timescale hence the maximal energy that can be reached; it 
is discussed in more detail in Sec. 13.11 For the time being, it suf- 
fices to note that the Larmor radius of the bulk of the electrons, with 
minimum Lorentz factor 7 m = \p— 1| _1 \p — 2|e e 7bim p /m e is so 
small that these electrons can only explore the undecayed part of 
the turbulence: 



r-L (7m) 



-V2 C ^ \ 
-o.3e B ,-2~ — ^ A P 



(9) 



The microturbulence thus controls the acceleration of the bulk of 
electrons, independently of how fast this microturbulence decays 
or how large the blast Lorentz factor may be. At Lorentz factors 
3> 7m of interest for high energy radiation, the electrons may start 
to explore the region fi(t) > 1. Then the transport becomes non 
trivial; its impact on acceleration is discussed in Sec. 13. II 

During the acceleration stage, the particle moves in a near 
ballistic manner and diffusive effects can be neglected, given that 
the return probability decreases fast with the number of steps of 
length Zscatt taken. Therefore, as a particle moves away from the 
shock on a distance scale x along the shock normal during a time 
t p ~ £'/(c cos 8), with 8 the angle to the shock normal, it explores a 
microturbulence that has decayed according to the laws given above 
with t ~ 3 cos 8 t p . This factor of 3 of course results from the con- 
vective velocity c/3 of the downstream plasma. 

On length scales much larger than i sca tt, a particle diffuses in 
the microturbulence and in a first approximation, one can describe 
its transport by advection with the downstream plasma. For such 
particles, t ~t p . The above slight difference between t and t p does 
not impact the results given further below and can be neglected in 
view of the uncertainties related to the time evolution of the turbu- 
lence. In the following, t and t p are thus be used interchangeably. 

The cooling history of an electron of Lorentz factor j e obeys 
the standard law 



dt 



1 SB(t) 2 (l + Y) 2 

— <TT 7e i 



()7T 



(10) 
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with Y the Compton cooling factor. For the time being, one con- 
siders the simple case Y -C 1; the influence of inverse Compton 
losses is discussed in detail in App. [Aland further below. One then 
defines a Lorentz factor 7 M + such that, if 7e > 7 M +, the particle 
cools in the undecayed microturbulence where n(t) < 1 (at time 
t < t M +), while if 7^ < 7 M +, the particle does not cool in that 
layer but further on. Writing the synchrotron cooling time for par- 
ticles of Lorentz factor "/ e in a (constant) magnetic field of strength 
SB as t syn [7 e ; SB], the Lorentz factor 7 M + can also be defined as 
the solution of 



up to a (comoving) hydrodynamical timescale td yn (see Eq, IA3b 



7m+ : M[*syn(7M+;£-B (J 



In terms of the fiducial values of interest here, 



7m+ 



3 x 10 9 4 /4 E 5 



1/4 -1/4 



(id 



(12) 



with A Mi 2 = A M /100. This value of 7 M + generally exceeds 
the maximal Lorentz factors that can be achieved through shock 
acceleration, therefore particles cool outside this undecayed tur- 
bulent layer, unless A M is larger than expected, as discussed in 
Sec. 13. II The latter may well happen, if for instance a\ > 2 and/or 
A„ » 10c/cj p i. 

If a particle exits the acceleration process with a Lorentz factor 
7e,o, then at time t, 



7e,0 



1 + 7s,o/7m+ 



7e,0 



1 + (1 + at)- 1 Ht) 1 + at + a t ] 7 £ ,o/7m+ 



(Mt) « i) 
(M*) » i) 

(13) 

If 7e0 > 7m+> me particle cools down to 7 M + within the layer 
where fj,(t) < 1 (i.e. t < i M +) and subsequently, it continues cool- 
ing if —1 < at < or stops its cooling if a* < — 1. If 7 e o < 7fl + 
on the contrary, the particle either cools later in the decaying mi- 
croturbulence if —1 < at < 0, or not if at < — 1. In any case, the 
particle of course eventually cools in the background shock com- 
pressed field (notwithstanding issues related to the available hydro- 
dynamical timescale). Whichever occurs influences the afterglow 
light curve and spectral energy distribution. 

If —1 < at < 0, it is convenient to define a second Lorentz 
factor, 7 M ~, as the Lorentz factor for which cooling occurs on a 
timescale t M _ such that SB = Bd, i.e. at the time at which 

the turbulence field has relaxed to the background shock com- 
pressed value Bd = 47bi-B u . For the time being, no considera- 
tion is made of the hydrodynamical timescale of the blast. Then, 
if 7m > 7/j- (and a t > —1), most particles cool in the decaying 
microturbulent layer. One finds 



7m- 
7m4 



-(l + "t)/«i 



(14) 



in terms of the upstream magnetization parameter a u = 
B^/(>iimmpC 2 ). The Lorentz factor 7 M _ depends exponentially 
on at; it may therefore take very different values. 

To calculate the radiative signature, one integrates over the 
cooling history of the electron population, as in Gruzinov & Wax- 
man (1999), although for simplicity, the calculation is done in a 
one-dimensional quasi-steady state approximation, meaning that 
the secular hydrodynamical evolution of the blast is neglected in the 
course of this integration over the blast width. This method corre- 
sponds to the steady state approximation of Sari et al. (1998), when 
calculating the stationary electron distribution in a homogeneous 
shell. The spectral power density of the blast can then be written in 
the downstream frame in terms of an integral over a particle history, 



P' = 



dN e d£, y „(7e.o) 

a d T e -° / dt a a+ > (15) 

d7e,o Jo dvdt 



with dE Byn /dudt the spectral power density radiated by an elec- 
tron at time t, of initial Lorentz factor 7 e ,o and of cooling history 
given by Eq. H3\ : see also Eq. lAll The above expression is folded 
over the injection distribution dN e /dy e ,o, which is assumed to take 
a power law form between 7 m and y max 2> 7m : 



4^^V P d 7e , 



with 



7m \ 7m 



13b + i J 4Trnr 2 c 



(16) 



(17) 



the number of electrons swept and shock accelerated by the shock 
wave per unit time, as measured in the downstream frame. 

The detailed calculation of the spectral power P v is carried 
out in Appendix [A] for different relevant cases. In particular, two 
distinctions have to be made: whether the microturbulence decays 
rapidly beyond i M + or not and whether inverse Compton losses 
contribute significantly to the cooling of electrons. These cases are 
examined in turn in the next subsections, in parallel to the discus- 
sion of App.lAl 



2.3 Gradual decay, no inverse Compton losses 

As discussed in App. [A] the gradual evolution of the microturbu- 
lence behind the shock affects the spectro-temporal flux F v in var- 
ious ways. For one, particles of different Lorentz factors cool in 
different magnetic fields, with particles of lower energy experienc- 
ing lower magnetic fields at cooling. This implies that the charac- 
teristic synchrotron frequencies are modified and more specifically, 
that ratios of characteristic frequencies are stretched with respect 
to the standard case of a homogeneous turbulent layer. Regarding 
the cooling frequency v c , the modification is non- trivial because 
the cooling Lorentz factor 7c itself depends in a non-trivial way on 
the temporal decay index of the turbulence, see Eq. lA14l The char- 
acteristic frequency u m associated to particles of Lorentz 7 m can 
also be modified in a non- trivial way, since i/ m = u p [ 7m ; SB lla ], 
with i/ p defined as the synchrotron peak frequency of particles of 
Lorentz factor 7 m in a magnetic field of strength <5£> 7m , see Eq. |A6l 
The magnetic field SB lm in which the particles of Lorentz factor 
7m radiate most of their energy takes the shock compressed value 
B d if both 7 m < 7m- and t dyn > but SB^t^/t^) 01 * /2 
if idyn < tp- and 7 nl < 7 c, or <5-B M (7m/7M+)~' 5t , with S t = 
Qt/[2(1 + at)] defined in Eq. lA17l The ratio tdyn/^M- determines 
whether the turbulence has relaxed to Bd by the back of the blast 
or not, and this value depends exponentially on at- 



tdyn 



1.6 x 10 4 e 6 ' 85/at S: 



2/a t 3/6+1/ a t 



xEi 



1/8 l/a t 



l t 5/S 



1 M.2 I '2 



-5/8 
-,0.3 



(18) 



In practice, it can take small or large values at different times, even 
for at = —0.5 for which the numerical prefactor becomes 0.018. 

In direct consequence of the above, the deceleration of the 
blast implies a non-trivial temporal evolution of the characteristic 
frequencies. This modifies the standard temporal evolution of F v . 
Furthermore, as a particle gets advected away from the shock with 
the microturbulence, it radiates at decreasing frequencies, whether 
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Figure 1. Comparison between two synchrotron spectra: in a homogeneous 
turbulence of strength eg = 10 -2 (top red curve), and in a decaying mi- 
croturbulence such that at = —0.5, A M = 10 2 . For both, the blast param- 
eters are 7b = 245, n_3 = 1 and the injection distribution index p = 2.2. 
The regime is slow cooling, and tdyn < *u— (corresponding to case 1 of 
Fig. lAll for the microturbulent model). 



it cools efficiently or not. Consequently, the changing magnetic 
field also modifies the spectral slope of the flux F v . 

Appendix [A] provides a detailed discussion of the possible 
cases, with detailed expressions for the characteristic frequencies 
Vm, v c , and t> M _ , depending on their respective orderings. Fol- 
lowing App. [A] the present discussion does not consider the cases 
in which either v m > or u c > Vn+, because such cases appear 
rather extreme in terms of the parameters characterizing the tur- 
bulence, as discussed in the former Section. Moreover, these cases 
tend to the standard model of a synchrotron afterglow in a homo- 
geneous turbulence when both ^ m > and v c > and they 
can be easily recovered from App. [A] 

Figure [TJpresents a concrete example of a spectral energy dis- 
tribution, comparing the standard prediction for a homogeneous 
turbulence with es = 1CP 2 (blue line) to a time evolving micro- 
turbulence with a± — —0.5, also starting at es = 1CP 2 (red line), 
with A M = 10 2 . Both models assume 7b = 245 corresponding to 
an observer time t b s — 100 s for a blast at z = 1 with £53 = 1, 
n_3 = 1, an injection slope p = 2.2 and a circumburst medium of 
constant density. Inverse Compton losses are neglected throughout 
the blast in this example. Figure[T]reveals the characteristic stretch 
of frequency range, with v m ~ 2 x 10 lf> Hz for the dynamical mi- 
croturbulent model (resp. u m ~ 2.3 x 10 17 Hz in the homogeneous 
turbulence) and u c ~ 1.1 x 10 22 Hz (resp. v c ~ 7.6 x 10 18 Hz). 
The frequency = 2 x 10 27 Hz lies outside the range of Fig.Q] 

For the above case of slow cooling, the spectral indices at 
low and intermediate frequencies, respectively p = 1/3 and /3 = 
-(1 - p)/2 (with F v oc t a u~ p ) remain unaffected in the pres- 
ence of decaying microturbulence. However the temporal index a 
is modified in these cases, even at low frequencies, which opens the 
possibility of testing such cases through the temporal behavior of an 
early follow-up in the optical. Section [01 below offers a compari- 
son of such spectra with the observed light curve of GRB090510. 
Regarding the fast cooling part of the electron population, both the 
spectral and temporal indices are modified, see Table [ATI for their 
detailed values. See also Fig. lAll for an illustration of possible syn- 
chrotron spectra. 

Synchrotron self-absorption is negligible at all frequencies 
shown in this figure, see Sec. lA5l for details. 




v[Hz] 

Figure 2. Same as Fig. [T] now comparing the synchrotron spectrum in a 
homogeneous turbulence of strength eg = 10 _ 2 (top red curve), with a 
decaying microturbulence such that at = —1.8, A M = 2.7 X 10 4 (bottom 
blue curve). For both, the blast parameters are 7b = 245, n_3 = 1 and the 
injection distribution index p = 2.2. The orange dashed line represents the 
secondary synchrotron component associated to cooling in the background 
shock compressed field (here B u = 10 fJ.G), which emerges on top of the 
microturbulent component because id yn = 1.3t M — (corresponding to case 
2 of Fig, lA2l . The regime is slow cooling. 



2.4 Rapidly decaying microturbulence (no inverse Compton 
losses) 

There are two main differences between the synchrotron spectra 
of gradually vs rapidly decaying microturbulence. In the former 
case, particles may cool in the decaying microturbulence layer, 
while in the latter, particles either cool in the undecayed region 
of short extent, if their Lorentz is sufficiently large, or in the back- 
ground magnetic field beyond the microturbulent layer otherwise, 
in the absence of inverse Compton losses that is. This implies in 
particular that there is no well defined cooling Lorentz factor for 
the microturbulence. Secondly, as the turbulence decays rapidly, 
most of the synchrotron power is emitted in the region of largest 
magnetic power, hence the flux vF v associated to the microturbu- 
lent layer peaks at At times such that tdyn > t^—, cooling 
in the background shock compressed field leads to the emergence 
of a secondary synchrotron component on top of the former, and 
one may now define a cooling Lorentz factor in the background 
compressed field. At frequency v m , the ratio between these two 
components is of order tdyn Am- > 1, at the benefit of the latter. 
At the exit of the microturbulent layer, the maximal Lorentz fac- 
tor cannot exceed 7 M +, so that the secondary component cuts-off at 
most at v v py M +; £>d], which falls short of the GeV range, see the 
discussion in Sec. 13.11 In practice, this suggests that most of the 
low energy emission in the optical and X-ray domains result from 
cooling in the background field, while the highest energy emission 
can be attributed to the presence of the microturbulence. 

Figure |2]presents a concrete example of a spectral energy dis- 
tribution, comparing the standard prediction for a homogeneous 
turbulence with es = 10 -2 (top red line) to a time evolving micro- 
turbulence with at — —1.8, also starting at e_e = 10~ 2 (bottom 
blue curve). The microturbulent layer is such that A M = 2.7 x 10 4 , 
corresponding for instance to a\ — 3, A M = 30c/oj p i. The up- 
stream magnetic field B u = 10 fiG and as in Fig. Q] n_3 = 1, 
7b = 245, t bs = 100 s, z = 1, which implies in particular that 
tdyn = 1.3t M _. Therefore, a secondary synchrotron associated to 
cooling in Bd emerges on top of the microturbulent component. 
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The corresponding characteristic frequencies for the microturbu- 
lent component are i/ m ~ 5.7x 10 15 Hz, u m ,SB tl — 2.3xl0 17 Hz 
and Vn-\- ~ 2.8 x 10 22 Hz; the cut-off frequency of the secondary 
synchrotron component v^o = 0.9 x 10 21 Hz. In such a scenario, 
the microturbulent component would dominate in the X-ray and 
at higher energies, while the secondary component dominates in 
the optical; at later times, the secondary component would come to 
dominate as well in the X-ray. 

Regarding the spectro-temporal evolution of the flux, one re- 
covers features similar to those discussed above in the case at > 
— 1. In particular, the slow cooling slopes 1/3 and (1— p)/2 remain 
unaffected, but the temporal index is different when t^ yYl < t M _, 
and for the fast cooling part of the electron population, both spec- 
tral and temporal indices are affected by at- The detailed values of 
a and j3 are given in Table lA2l in Appendix [A] See also Figs. IA2| 
and lA3l for an illustration of the spectral shapes of the synchrotron 
spectra. 

Synchrotron self-absorption is negligible at all frequencies 
shown in this figure, see Sec. lA5l for details. 



2.5 Strong inverse Compton losses 

Accounting for inverse Compton losses modifies of course the cool- 
ing history of the particle. The importance of inverse Compton 
losses is generally quantified through the Y Compton parameter, 
which may be written in a first approximation (e.g. Sari & Esin 
2001, Wang et al. 2010) as 



Y(i) [1 + Y( 7 )] 



vFy [^kn(7)] 



./"cool 



(19) 



€ B vF u (z/ pcak ) 

where the notation Y (7) indicates that Y in general depends on the 
Lorentz factor of the particle, due to Klein-Nishina effects. These 
latter are quantified through the ratio of vF v [^kn(7)], i.e. the syn- 
chrotron flux at the frequency at which Klein-Nishina suppression 
becomes effective, to the peak of the synchrotron flux, vF v (^ pea k). 
Finally, the factor / 0oo i denotes the fraction of cooling electrons, 
~ 1 for fast cooling, ~ (jc/jm) ~ s for slow cooling. The double 
dependence of Y on Lorentz factor 7 and distance to the shock 
front (through SB) that arises in the present model renders this 
problem quite complex. To simplify this task, it is assumed here that 
inverse Compton losses dominate everywhere throughout the blast. 
This would occur for instance, if e e > es in the undecayed part 
of the microturbulent layer (a generic assumption), if / coo i does 
not lie too far below unity, and at energies such that Klein-Nishina 
effects can be neglected. Under such conditions, one may expect 
Y > 1 everywhere else, because the magnetic field decays away 
from the shock. 

In the presence of inverse Compton losses, one can always de- 
fine a cooling Lorentz factor 7 C irrespectively of the value of at- 
As the cooling history of the particles is modified, so are the spec- 
tral slopes in the fast cooling part if at > — 4/(p + 1), and also 
at low frequencies if at < — 4/ (p + 1). In this latter situation in- 
deed, the synchrotron spectrum is dominated at high frequencies by 
a slow cooling spectrum in the undecayed microturbulence, with 
a low energy extension from f m to f m ^s M [see Eq. IA30I with 
slope 1 + 2 /at. At low frequencies, a secondary component as- 
sociated to cooling in the background compressed field emerges 
if both tdyn < ty,- and 7 m < 7 M -. The synchrotron spectra for 
at > — 1 share features similar to those discussed in Sec. 12.31 up 
to the modifications of the spectro-temporal indices. One notewor- 
thy distinction is the fact that the fast cooling spectral index has 
become — (p + at /2)/(2 — at /2), which may become larger than 
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Figure 3. Same as Fig. ff] now comparing the synchrotron spectrum in a 
homogeneous turbulence of strength es = 10 (top red curve) including 
inverse Compton losses with Y = 3, with a decaying microturbulence such 
that at = —0.8, A^ = 10 2 , also including inverse Compton losses with 
Y/j, = 3 (bottom blue curve). For both, the blast parameters are 71, = 245, 
n_3 = 1 and the injection distribution index p = 2.2. Here id yn <C t^— 
(corresponding to case 1 of Fig. lA4l . The regime is slow cooling. 



— 1 if at < 2 — p, in which case most of the power lies at . This 
may be attributed to a scaling of the Compton parameter with fre- 
quency: higher frequencies correspond to higher Lorentz factors, 
hence to cooling in a region of higher magnetic field, where the 
Compton parameter is smaller, which implies that a smaller frac- 
tion of the particles energy is channeled into the inverse Compton 
component. 

Figure |3]presents a concrete example of a spectral energy dis- 
tribution, comparing the standard prediction for a homogeneous 
turbulence with £b = 10 -2 including inverse Compton losses 
with Y — 3 (top red line) to a time evolving microturbulence 
with at = —0.8, also starting at es = 10 -2 , with A M = 10 2 
and y M = 3 (bottom blue curve). For this case, idyn <C tp-, 
corresponding to case 1 of Fig. IA4I other parameters remain un- 
changed compared to previous figures. The corresponding charac- 
teristic frequencies for the microturbulent component are v m ~ 
4.7 x 10 15 Hz, v c ~ 1.0 x 10 16 Hz and k„ + ~ 1.3 x 10 26 Hz. 

The spectro-temporal indices a and /3 are given in Table lA3l 
forat > —4/(0+1) and lA4l in the opposite limit at < — 4/(p+l). 
The generic spectral shapes are illustrated in Figs. I A4 1 and I A 5 1 for 
these two cases, respectively. 

Here as well, synchrotron self-absorption is negligible at all 
frequencies shown in this figure, see Sec. |A5| for details. 



3 DISCUSSION 

Most of the discussion so far has ignored the characteristic fre- 
quency associated to the maximal energy of the acceleration pro- 
cess. Yet, the production of GeV photons through synchrotron ra- 
diation of electrons already push ultra-relativistic Fermi accelera- 
tion to its limits, as discussed e.g. in Piran & Nakar (2010), Kirk & 
Reville (2010), Barniol-Duran & Kumar (2010, 2011a), Lemoine 
& Pelletier (2011c), Sagi & Nakar (2012), Bykov et al. (2012). 
These studies assume a homogeneous downstream turbulence, with 
either microscale high power turbulence or simply a shock com- 
pressed magnetic field. Section [3~T| extends these calculations to 
an evolving microturbulence, which brings in further constraints 
and new phenomena. A brief comparison with the light curve of 
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GRB090510, which so far provides the earliest follow-up in the 
optical through GeV, is then proposed. 

3.1 Acceleration to high energies 

The maximal energy is determined through the comparison of the 
acceleration timescale t acc to other relevant timescales, e.g. the en- 
ergy loss and dynamical timescales. Given that the energy gain per 
Fermi cycle is of order 2 (Gallant & Achterberg 1999, Achterberg 
et al. 2001, Lemoine & Pelletier 2003), it suffices to compare in 
each respective rest frame the downstream £ rcs | d and the upstream 
t r cs | u residence times to the other timescales. 

3.1.1 Upstream 

The upstream residence time t ros | u ~ 5tL,o|u/7sh (Lemoine & 
Pelletier 2003) - assuming that the particle only interacts with 
the background field - where <l,o = 7b7eWi e c/ (eB u ) and 7 e 
represents the particle Lorentz factor in the downstream frame, 
7sh = v^lb- Considering this residence time leads to a conser- 
vative limit on the maximal energy, because scattering in the mi- 
croturbulent field seeded in the shock precursor should also play a 
role, see Plotnikov et al. (2012) for further discussion on this issue. 
Comparing the above residence time to the age of the shock wave 
r/c then leads to a maximal Lorentz factor (downstream frame) 

z> 1.2 x 10 9 fl_ B E& 4 n:S /4 # 4 *T#3 . (20) 

Depending on the dynamics of the microturbulence downstream, 
this maximal Lorentz leads to to photons of energy 

t 7,max — 1 ieV-D_5£/ 53 fl_g t 2 Z +,0.3 £ B,-2> ^- L > 

through radiation in the 8B^ field, or 

e^Lx * 30GeVB 3 _ 5 E 3 5 i 4 nZl /4 h 1/4 z-%i , (22) 

if the particle radiates in the background shock compressed field 
Bd, as proposed in Barniol-Duran & Kumar (2010, 2011a). This 
latter case applies in particular if a t < — 1, Tmax < 7n+ and if 
downstream inverse Compton losses can be neglected. If however, 
7m2x > 7m+> then the particle radiates in the stronger <5S M , while 
if a t > —1 and 7m2 x < 7n+> the particle radiates in a field of 
strength SB ~ SB^ ^7max/7/j+) . which for these high ener- 
gies lies in practice close to SB^. For instance, if at = —0.5, one 
finds 

( a ) ^ nKTirD 2 ^!? -1/4,-0.5 -0.5 a 0.5 

(23) 

Equation[22]matches the "confinement" estimates of Piran & Nakar 
(2010), Barniol-Duran & Kumar (2011a), who considered only a 
shock compressed background field in the downstream, while the 
other estimates Eqs.|2T]|23]hold in the present more realistic setting 
because of the high value of magnetic field behind the shock front. 

Note that the particle is actually confined close to the shock 
front, because in the ultra-relativistic regime, the shock front never 
lags further than ct res i u /7gh behind the particle; escape can only 
take place through the boundaries, but this would lead to a less 
stringent condition on maximal energy than the above age con- 
straint, see e.g. Bykov et al. (2012). 

The above suggests that the upstream magnetic field can lie 
well below 10 /iG and yet lead to the production of multi-GeV pho- 
tons because the maximal Lorentz factor particles cool in a strong 



field, close to the microturbulent value SB^ behind the shock. Fur- 
thermore, if one neglects altogether the upstream magnetic field, 
the discussion in Plotnikov et al. (2012) suggests that scattering the 
microturbulent field generated in the shock precursor may lead to 
Lorentz factors sufficiently large to produce GeV photons through 
synchrotron radiation in SBfj,. 

Li & Waxman (2006) and Li & Zhao (201 1) have argued that 
the comparison of the upstream residence time with the timescale 
for inverse Compton losses imposes a stringent lower bound on 
the upstream magnetic field, suggesting an efficient amplifica- 
tion of the latter. More specifically, Li & Zhao (2011) find that 
1 mG n 9 /s < B u < 100 mG n /s (with B u denoting the total 
magnetic field upstream of the blast) is imposed by the synchrotron 
model of Fermi— LAT bursts with extended GeV emission. These 
bounds have been challenged in Bamiol-Duran & Kumar (2011b) 
and Sagi & Nakar (2012). Nevertheless, such bounds agree well 
with the general picture that microinstabilities self-generate a mi- 
croturbulent field in the shock precursor, which is then transmitted 
downstream and there forms the microturbulent layer. One expects 
that in the shock precursor, the microturbulent field B u ~ SB^/^ 
(as measured in the upstream frame), in which case one finds 
B u ~ lmGegij"-'?- This field satisfies the constraints of Li 
& Zhao (2011). ' 

3.1.2 Downstream 

Consider now the situation downstream, which leads to much more 
severe contraints on e 7 , max through the comparison of £ rGS | d with 
the timescale of energy loss. Out of commodity and simplicity, the 
limit on the maximal energy is often quoted with a Bohm estimate 
t reB |d — kb*l, with = (j e m £ c) / (eSB) in the downstream 
frame. Ignoring inverse Compton losses, this leads to 

- 2.1 x 10\-^r 6 E; 3 1/16 n-_ 3 3 /16 e-%t , (24) 

assuming for the moment that SB is homogeneous, with strength 
characterized by e_e • This leads to 

4 6) max =a 30GeV K B 1 £5 3 / V / VJ /8 , (25) 

which matches the estimate of Piran & Nakar (2010) up to a factor 
2. However, this estimate must be corrected for two effects: (1) the 
Bohm approximation likely breaks down at these large Lorentz fac- 
tors; (2) the evolution of the microturbulence away from the shock 
front modifies the scattering rate of the particles. 

The Bohm approximation for t ros | d fails at the maximal 
Lorentz factors because 

^ J ~ 1.7 x 10 2 X^tT 6 E^n-_\' le etL% » 1 , 
-V 

(26) 

for a Lorentz factor given by Eq. [24] This implies that the parti- 
cle only suffers a random small angle deflection of order A^/rL 
as it crosses a coherence cell of size A M = 10A Mi ic/oj p i, hence 
the residence timescale is rather given as: £ rcs | d ~ k sc (fl/A^) tl, 
much larger than the above Bohm estimate. This of course reduces 
significantly the maximal energy of photons (Kirk & Reville 2010, 
Lemoine & Pelletier 2011c, Plotnikov et al. 2012). At the present 
time, one cannot exclude that the spectrum of turbulent modes ex- 
tends to wavenumbers k ~ r^ 1 that would allow gyroresonant in- 
teractions in the downstream; whether or not this happens depends 
on the physics of instabilities in the far upstream, which are still 
debated. Were this the case, the Bohm estimate should nevertheless 
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be corrected by a factor that accounts for the diminished magnetic 
power at those resonant scales, compared to the maximum power 
at A^ 1 . In order to obtain a conservative estimate of the maximal 
energy, one should therefore rely on the above £ rcs | d oc r£. Com- 
paring this residence timescale with the synchrotron loss time, this 
would lead to a maximal Lorentz factor 



~(<0 

/max 



(27) 



3.9 x 10 7 k~ 1/3 Abulia 76 , 
corresponding to a maximum photon energy 

41 ax * 0.9GeV^ 3 EU\zl /12 Xl{!^ 2 _ 2 t-^ , (28) 

provided the radiation takes place in the microturbulent field of 
strength SB^. If radiation were to take place in the background 
shock compressed field Bd, the maximal photon energy would 
rather be 



30 MeV k: 



2/3 EU 4 n- 7/12 ^ 2/:i 



-3/4 



(29) 



However, once the evolution of the microturbulence is ac- 
counted for, these estimates are modified towards more optimistic 
values. Consider first the scaling of rx/A, with ri, now design- 
ing the gyroradius in the local magnetic field, and A the coherence 
length of this local magnetic field. Then 



— oc t 
A 



a t /2-l/a\ 



(30) 



With a\ ~ 2 — 3, the following qualitative picture emerges. If 
the turbulence decays gradually, meaning a t > — 1, the coherence 
length increases more rapidly away from the shock front than the 
magnetic field loses power, so that the scattering timescale tends 
towards a Bohm scaling (albeit, in a weaker magnetic field) as time 
increases (see also Katz et al. 2007 for a similar picture). One then 
obtains (see further below) maximal energy estimates close to the 
Bohm scaling. If, however, the turbulence decays fast, meaning 
at <■ — 1, then the scattering timescale increases with time, pos- 
sibly faster than t if a t < — 1 — l/a\, in which case scattering be- 
comes impossible in the decaying part of the microturbulent layer. 
Furthermore, mirror like interactions in the background shock com- 
pressed magnetic field, which lies transverse to the shock normal, 
do not allow repeated returns to the shock front (Lemoine et al. 
2006, see also Sec. |3.3t . The particles must then scatter in the un- 
decayed microturbulent field, if at < — 1, in which case the maxi- 
mal energy is then determined by the comparison of t rcs |d to 
These two scenarios are examined in turn. 

Assume first at > — 1, meaning more specifically at > 



-2/ ax- Then the small angle scattering timescale t 



t 



-at — 



rcs|d 



which must be compared to the synchrotron loss time 
t syn oc 7mi x 5B M {t le s\d/tfj,+) where the last factor accounts 
for the evolution of the magnetic field strength with time, outside 
the undecayed layer. This leads to 

\2-\ -(l + a t )/(3+3Q t + l/a A ) 

(31) 



/max 



7^+ 



A,, ct 



The term within the brackets can be rewritten in a more compact 
way as 47^ + r e A^ 1 /9, with r e the classical electron radius, so that 
in the limit ax — > +oo, one recovers the estimate 7 



mogeneous microturbulence. The fiducial values at = 
ax — 2 discussed previously lead however to 



( c ) • U 

Jn,i x in a ho- 
0.5 and 



/max 



1.2 x 10 



a -0.25 t-,-0. 06 - 
A M,2 #53 n - 



0.19 - 

3 e, 



0.25 -,0.25,0.19 

,-2 A M,1 r 2 , (.•>*■) 



which, through synchrotron radiation in a magnetic field of strength 
SB ~ SB^ ^7m2x/7/j+ ) leads to photons of typical energy 



2GeV£ 5 °3 22 A2' 6 1 3 A; 



0.094 0.38 .- 
3 e B,-2% 



0.(i(i 



(33) 



One should stress that this result depends on the particular value of 
a t and ax, but for the above fiducial values, it remains close to the 
Bohm estimate and allows the production of GeV photons for con- 
servative assumptions. Note finally that the stretching of the coher- 
ence length with time does not imply that the particles of Lorentz 
factor 7m2x suffer gyroresonant interactions away from the shock, 
nor does it require that the turbulent power spectrum extends over 
many decades. For the above fiducial values, the particles interact 
with modes of wavelength / ~ 30 A M -C 7%. 

In a fast decaying scenario, meaning more exactly a t < 
— 1 — 1/ax, scattering must take place in the undecayed part of 
the microturbulent layer, otherwise the particles would not see the 
microturbulent layer and Fermi acceleration would not take place. 
The comparison £ rcs | d < t M + leads to 



-Y (e) 

/ max 



4xlO e .- 1/2 ^ / X{i 2 4 / > 
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which would lead to O(10) MeV photons only, even through radi- 
ation in a field as strong as SB^. To reach the GeV range, ceteris 
paribus, one needs to increase the spatial extent of the undecayed 
microturbulence, i.e. to increase A Mi 2 by 2 — 3 orders of magni- 
tude. For A M > 10 4 , indeed, 7max > 7m+. so that the particle both 
scatters and cools in the undecayed part of the microturbulent layer. 

To summarize, the above analysis of scattering and maximal 
photon energies leads to the following qualitative picture. If the 
turbulence decays gradually, the particles scatter while the micro- 
turbulence decays and the maximal energy of the photons typically 
falls in the GeV range for fiducial parameters and rather conser- 
vative estimates of the scattering properties. If, however, the tur- 
bulence decays rapidly, GeV photons can only be produced if the 
undecayed part of the turbulent layer extends far enough to acco- 
modate both the scattering and the cooling of the maximal energy 
electrons in that layer. This latter requirement places strong con- 
straint on the parameter A M that characterizes for how long the un- 
decayed turbulence can survive. Requisite values can be obtained 
if, for instance, A M ~ 30c/oj p i and ax ~ 3, which cannot be ex- 
cluded at present. Finally, it has also been noted that the constraints 
on the upstream residence time are much weaker if the maximal 
energy electrons can radiate in a field close to 5£? M . In practice, 
this implies that the value of the upstream magnetic field is not 
well constrained. This alleviates the apparent strong magnetization 
problem discussed in He et al. (201 1). 



3.2 Early GRB light curves 

Section [2] together with Appendix [A] have emphasized the differ- 
ences between the standard synchrotron spectra for homogeneous 
turbulence and those calculated with the account of decay of the 
microturbulence. This Section confronts such signatures with the 
observational data on GRB090510, which so far has provided the 
earliest follow-up on a broad spectral range, with near simultane- 
ous detection in the optical, X range and GeV range (de Pasquale 
et al. 2010; Ukwatta et al. 2009; Guiriec et al. 2009; Ackermann 
et al. 2010; Longo et al. 2009; Hoversten et al. 2009; Golenetskii 
et al. 2009; Ohmori et al. 2009). This choice is further motivated 
by the analyses of Barniol-Duran & Kumar (2009, 2010, 2011a) 
and He et al. (201 1), who have argued that the extended emission, 
from the optical to the GeV range is compatible with an afterglow 
spectrum calculated without magnetic field generation beyond di- 
rect shock compression of the interstellar field. One may also in- 
terpret the data with a high magnetization of the blast, closer to the 
traditional estimates of es, although it then requires an extraordi- 
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narily low upstream density < 10 6 cm 3 (e.g. de Pasquale et al. 
2009, Corsi et al. 2010); this possibility is not considered here. The 
model of Barniol-Duran & Kumar (2009, 2010, 2011a) then sug- 
gests that magnetic field generation processes are ineffective, and 
that this burst provides the closest relative to a "clean" relativistic 
blast wave in a weakly magnetized medium 

As discussed in the previous Sections, one should expect that 
the GeV photons have been produced in the microturbulent layer 
rather than in the background magnetic field. However, lower en- 
ergy photons in the optical and in the X range are typically pro- 
duced by particles of Lorentz 7 m , which lose energy on a much 
longer timescale than the particles of Lorentz factor 7 max , hence 
in a weaker magnetic field, possibly the background shock com- 
pressed value Bd- In any case the microturbulence is to play a key 
role in shaping the spectra from the GeV down to lower energies. 
This motivates further the search for a signature of this microturbu- 
lence. 

Unfortunately, the comparison is not straightforward, because 
of the diversity and complexity of the synchrotron spectra. The 
quantities that characterize the microturbulence appear as new free 
parameters from the point of view of phenomenology, hence they 
enlarge the dimensionality of parameter space and add degener- 
acy. For this reason, the following does not attempt to fit accu- 
rately the light curves in the GeV, X-ray and optical, but rather 
accomodates the different spectro-temporal indices in the differ- 
ent time domains of the observations. The key observation of the 
Barniol-Duran & Kumar (2009, 2010, 2011a) interpretation is that 
the spectrum is produced in a slow cooling regime, and that in the 
case of GRB0905 10, v m transits across the optical at 10 3 s, in order 
to produce the observed break. Once this condition is satisfied, and 
provided spectro-temporal indices match the observed values, one 
obtains a satisfactory fit to the light curves (as it has been checked). 
The following discussion reveals that the current data do not al- 
low one to pick out a preferred model, meaning conversely that it 
is possible to reproduce the salient features of the Barniol-Duran 
& Kumar (2009, 2010, 2011a) model with concrete characteriza- 
tions of the microturbulence. Each model leads to rather specific 
signatures that might be identified with further higher accuracy and 
earlier observations. Some of these signatures are pointed out in the 
following, although more work is needed to define an observation 
strategy capable of pinpointing the characteristics of the microtur- 
bulence. 

GRB090510 is a short burst at redshift z = 0.92 whose light 
curve is characterized as follows: emission is seen in the LAT 
range up to ~ 100 s, with spectro-temporal indices qgcV — 1-4, 
/3ggV — 1-1; X-ray and optical follow-up start at about this time; 
from ~ 100 s to about 1.5 x 10 3 s, q x ^ 0.74, increasing to 2.2 
afterwards, while /3x — 0.5 — 0.8; a opt — —0.5 before 10 3 s, 
then 1.1 afterwards. The late time evolution beyond 10 3 s is more 
difficult to reproduce, in particular the steepening of the X-ray 
lightcurve, which is attributed by He et al. (2011) to sideways ex- 
pansion, but which does not match the apparent shallower decay in 
the optical up to 2 x 10 4 s (Nicuesa Guelbenzu et al. 2012). 

To compare this light curve with the afterglow computed in 
App. [A] consider first the possiblity of a t > —1 (gradual decay) 
without inverse Compton losses. The latter assumption has been 
discussed at length in He et al. (2011). As shown in Fig. lAll the 
spectral shape of the synchrotron spectrum in the slow cooling limit 
with at > — 1 retains the same form as the standard afterglow, but 
the values of the characteristic frequencies and therefore the time 
dependences are modified. Therefore, /3x — 0.6 if ^ m < 1 keV < 
v Q as observed. In the model of Barniol-Duran & Kumar (2010), 



He et al. (201 1), u c lies below the LAT range, so that 

p- at (l-2p)/2 
PaeV ~ 2 + 3a t /2 ' (35) 

which takes values between —1.10 and —1.00 for p = 2.2 and 
at £ [0, —1], in good agreement with the data, although the in- 
fluence of at is too weak to provide significant constraints. The 
temporal slope in the GeV range also depends mildly on (Xf, tak- 
ing values from a = 1.15 at at = (p = 2.2) to a — 1.25 
at at = — 1 (p = 2.2). It fits with the observations, given that the 
measured slope may be contaminated by prompt GeV emission, see 
the discussion in He et al. (2011). The best probes for a t remain 
the temporal slopes in the X and optical domains: a takes values 
between 0.9 and 1.4 in the X range (for p = 2.2), and between 
-0.50 and -0.30 (p = 2.2) in the optical when 1 eV < v m . Cur- 
rent data do not allow a precise determination of at, but it may be 
noted that the fiducial value a t = —0.5 fits well the data. Finally, 
requiring v m ~ 1 eV at 10 s determines a t as a function of the 
other parameters. For the fiducial values characterizing the turbu- 
lence, n ~ 10 -3 cm -3 and B d < 1 /iG, v m = u p [y m ; SB(t dyn )] 
because the turbulence has not had time yet to relax to the back- 
ground magnetic field B d . Then 
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so that v m ~ 1 eV at 10 s implies at = —0.6 up to small loga- 
rithmic corrections. It is somewhat remarkable that the value of at 
falls close to the fiducial value discussed previously. Also, one may 
note that the upstream magnetic field does not enter Eq. [36] be- 
cause the particles actually cool in the decaying microturbulence. 
One may verify that the above choice of parameters gives a light 
curve in good agreement with the observed data, except for the late 
time shallow decay in the optical mentioned above. 

To pursue this comparison, consider now the possibility that 
— 4/(p +1) < at < —1. The fast decay of the microturbulence 
in this case actually mimics somewhat the scenario originally pro- 
posed by Barniol-Duran & Kumar (2010, 2011a): the particles get 
accelerated in a microturbulent layer behind the shock front but 
cool where the microturbulence has died away. Although, as dis- 
cussed in some details in Sec. 13. II it appears necessary to require 
that the particles of Lorentz factor 7 max actually cool in the micro- 
turbulent layer in order to produce the GeV photons, in other words 
7m+ ~ 4 x 10 7 at tobe ~ 100 s, which implies 

A ^ ln 4 c ,-l/4 -1/4 -1 ,1/4 

or to put it more simply, that the undecayed part of the microturbu- 
lent layer with 63,-2 ~ 10 -2 extends for some 10 4 skin depths at 
least. 

Two different behaviors can be observed, depending on the ra- 
tio tdyn/ip-, i-e. whether the turbulence has relaxed down to B c \, 
in which case a synchrotron component associated to Bd would 
emerge on top of the microturbulent synchrotron spectrum. The 
maximum value of idyn/ip- is obtained for a t — > — 4/(p + 1) ~ 
-1.3: 
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so that both possibilities have to be envisaged. If t c \ yn < t^- up 
to late times, and i> m < 1 eV < 17,0 beyond 10 3 s, the optical 
flux should decay fast with a op t — 1-4 — > 1.6, which disfavors 
this possibility. If now tdyn > t^- at some point, the synchrotron 
component produced by the cooling of particles in Bd emerges and 
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dominates over that produced by the microturbulence, at least at 
low frequencies (optical, X). This Bd component cuts off at v^q, = 
v-p [7n+;Sd], as discussed in Sec. 12.41 f M o increases with time, 
and for reasonable choices of parameters, it lies above the X range. 
At the time at which td yn = t^-, however, the temporal decay 
slope in the X range changes from (6p — l)/8 (due to cooling in 
the microturbulence) to 3(p — l)/4, which implies a change from 
steep to shallow, which is not observed. Therefore, this scenario is 
disfavored as well. It is difficult at the present to rule it out clearly, 
as other effects, due to jet breaking for instance, could complicate 
the temporal scaling. 

Consider now the possibility —3 < at < — 4/(p + 1), with- 
out inverse Compton losses as above. Similar considerations re- 
garding the maximal energy suggests that 7 M + < 4 x 10 7 , imply- 
ing A M > 10 4 . The turbulence decays fast so that, for a reason- 
able choice of parameters, e.g. at < —1.5, ri-3 ~ 1, B-g ~ 1, 
the condition td yn > t^- may be fulfilled at t b s > 100 s. This 
implies that the background synchrotron component emerges on 
top of the microturbulent synchrotron spectrum, and dominates at 
low frequencies. One then recovers the standard afterglow spectrum 
proposed by Barniol-Duran & Kumar (2010), except that the GeV 
photons are produced in the microturbulent layer. In the GeV range, 
the indices (a, /3) match the standard predictions for the fast cool- 
ing population and thus agree with the data; the overall afterglow 
provides a satisfactory fit to the available data. In order to discrim- 
inate this model, one would need to access the early light curve in 
the frequency range v m < v < (i.e. X-ray) in which both 

the spectral and temporal indices evolves strongly with at'. (a,/3) 
go from (1.6, 0.5) to (0.1, —0.4) as at goes from —1.3 to —3. For 
reference, Vm,6B^ — 2.3 x 10 17 Hz E^eg 2 _ 2 e l-o.3t2~ 3 '' 2 - 

Another interesting solution arises in this scenario when 
tdyn < tft- up to late times. This may be realized if B_5 <C 1 
for instance. Then the frequency v m — v p [y m ; SB (tdyn)] transits 
through the optical at 10 3 s if at ~ —1.5, as required. During the 
interval 10 3 — 10 4 s, the optical lies in the range v m — ^ m ,is„ with 
(a, /3) = - [(24 + 7a t )/(Sa t ), 1 + 2/a t ]. To reconcile a with 
a op t ~ 0.8 (Nicuesa Gulbenzu et al. 2012), one would need to 
impose a t ~ —1.8, while /3 op t ~ 0.9 at that time rather requires 
a t — —1.1. Thus a t — —1.5 provides a compromise, which is 
marginally acceptable by the data, that presents the clear advan- 
tage of a shallower decay in the optical than in the X-ray at a time 
at which it is observed. At 10 4 s, v m ,sBu transits across the optical 
(see its expression above), so that the optical and the X-ray domains 
decline alongside at later times. 

Whichever ratio tdyn/t^- is considered, the model with a t < 
— 4/(p+ 1) predicts a clear signature in the early optical lightcurve, 
with Qf op t = —1/8 when 1 eV < v m , which leads to a mild 
growth the optical flux. De Pasquale et al. (2010) report a opt = 
— 0.5lg'jJ between 10 2 s and 10 3 s, while Nicuesa Gulbenzu et al. 
(2012) rather indicate a opt = —0.2 ± 0.2, therefore it is not possi- 
ble to rule out or confirm this model at present. 

Consider finally the possiblity of strong inverse Compton 
losses. As before, it is assumed here that these losses dominate 
throughout the blast, although a more realistic model should in- 
clude the possibility of weak inverse Compton losses at the highest 
energies due to the Klein-Nishina suppression of the cross section. 
Nevertheless, one may constrain this scenario with the optical and 
X data, as follows. 

The dominance of inverse Compton losses for particles of 



Lorentz factor y m implies that 7 C ~ 7 m , since 

7 C ~ i r r -l/2 -1/2 -1 -1 4 o.5 , Qm 
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see Eq. lA35l for the definition of Y^, Consider first the possibility 
at > — 4/(p + 1), this implies that i/ m < 1 keV and v c < 1 keV 
at tobs > 100 s, so the X range always fall in the fast cooling 
portion. The spectral slope cannot be — p/2 in that range, other- 
wise it could not match the observed /3x\ the slope must thus be 
-(p + a t /2)/(2 - a t /2), which fits fix ^ 0.5 - 0.8 only for 
at ~ —0.8 — ¥ —1.3. In turn, this would imply a positive value of 
a op t if v m > Vc (fast cooling), so the regime at time t b s > 100 s 
must be that of slow cooling (case 1 or 2 depicted in Fig.|A4j. One 
should further require that v m goes into the optical at 10 3 s as be- 
fore. These features can be achieved for the typical fiducial values 
considered before, n_3 ~ 1, B53 ~ 1, with a t ~ —0.8 — > —1.3. 
An interesting feature of this scenario is that the ratio 

^dyn / 1 

ix— can take values smaller or larger than unity for stan- 
dard parameters of the microturbulence and at ~ —1. As it de- 
pends sensitively on the value of the upstream magnetic field, as 
B u 2 ^ at , the transition from tdyn < t/i- (case 1 in Fig. IA4t to 
tdyn > tfi- (case 2 in Fig. lA4t may well occur around 10 2 — 10 3 s. 
Once tdyn/t,j- > 1, the optical domain lies in the fast cool- 
ing region corresponding to v c < 1 eV < 17,- and indices 
(a,/3) = [(3p-2)/4,p/2] ~ (1.1, 1.1). The optical domain 
then decays more slowly than the X domain, a being in marginal 
agreement with the inferred value a opt = 0.8 ± 0.1 (Nicuesa Gul- 
benzu et al. 2012). Note that the frequency i/ M _ increases in time as 
t 3 { 4 , but remains below the X domain if B-5 < 1. Therefore the 
break at 2 x 10 4 s in the optical cannot be attributed to another fre- 
quency crossing the optical domain; it might however result from 
jet sideways expansion in that model. In order to probe and dis- 
criminate this model, one would need to access the very early light 
curve, while it is in the fast cooling regime and measure the tem- 
poral decay slopes of the X and optical domains. The above mod- 
elling has implicitly assumed a Compton parameter Y M ~ 3; one 
can check that for such a value, the inverse Compton component of 
optical photons at 100 s remains below the synchrotron flux in the 
GeV domain. 

If at < — 4/(p + 1), the spectra depicted in Fig. \A5\ gen- 
erally resemble at high energies a slow cooling spectrum in the 
microturbulent field. Since i^ m .^B M iS 1 keV at t G b s > 100 s, the 
X domain always lies on the slow cooling portion of the microtur- 
bulent synchrotron spectrum with therefore correct values for ax, 
/3x- Among the 5 cases shown in Fig. |A5l the only two reasonable 
ones are 5 and 2. Indeed, case 1 is unlikely because it requires a 
small external density to achieve slow cooling, but then tdyn/tp- 
becomes large and case 2 actually applies. Similarly, case 4 is un- 
likely because v m ~ v c as argued above. Case 3 leads to a opt > 0, 
which goes contrary to the observational data. Therefore, cases 5 
and/or 2 remain as the viable alternatives. However, one must now 
impose !/ M _ < 1 keV at 10 2 — 10 4 s, so that the secondary syn- 
chrotron component associated to Bd cuts off below the X range, 
otherwise the X domain would fall in the fast cooling part of that 
synchrotron component. This constraint is fulfilled for the previous 
fiducial values B-5 < 1, n_3 ~ 1, E53 ~ 1 and for at > —1.5. 
Interestingly, the optical falls on the fast cooling part of the sec- 
ondary synchrotron component at > 10 3 s, since v m ~ v c , there- 
fore it decays less fast than the X range, with a in marginal agree- 
ment with the observed value, as previously. The temporal slope 
in the optical at very early times < 10 2 s provides an interesting 
check of this scenario: it should be flat or even weakly decaying, 
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corresponding to case 3 in Fig, IA5l which represents the analog of 
case 5 when tdyn < ip-. Finally, as discussed before, one must 
require A M > 10 4 in order to produce GeV photons directly in the 
undecayed part of the microturbulent layer. In this model, the mi- 
croturbulent layer produces both the GeV and X photons, while the 
optical photons are produced in the shock compressed background 
field. 



3.3 Further considerations 

According to the modelling of Barniol-Duran & Kumar (2009, 
2010, 201 la), the magnetization of the blast wave can take a large 
range of values, from eb ~ 10~ 9 to es ~ 10~ 3 , depending on 
the external density. The latter value is given in de Pasquale et al. 
(2009) and Corsi et al. (2010), assuming very low upstream den- 
sities n < 10~ 6 cm~ 3 . The general trend is that, independently 
of n, the downstream magnetic field appears to coincide with the 
shock compressed magnetic field = 47b-B u , corresponding to 
S u ~ 1 - 10 /iG. As noted in He et al. (2011), this generally 
implies large upstream magnetizations, since a u ~ es/2 if e_B is 
calculated in terms of £>d, whereas the typical interstellar magne- 
tization is more of the order of 10 -9 . As discussed in the previous 
sections, accounting for a decaying microturbulent layer alleviates 
the need for a large value of B u , which allows to reconcile the ap- 
parent e b with a more standard upstream magnetization. 

Nevertheless, it is interesting to ask what would happen if 
the upstream magnetization were truly high, e.g. if B u > 1 /iG 
and n <C 10~ 3 cm -3 . If a u is brought upwards by several orders 
of magnitude, then magnetization effects are to affect the shock 
physics. For instance, the PIC simulations of Sironi & Spitkovsky 
(2011) indicate that, for 7b = 15 and a u > 10~ 4 , the shock is no 
longer mediated by the filamentation instability, but by the mag- 
netic barrier associated to the (transverse) background magnetic 
field. As discussed in Lemoine & Pelletier (2010, 2011a), such a 
transition occurs when the growth timescale of the filamentation 
instability becomes larger than the timescale on which upstream 
plasma elements cross the shock precursor. This transition is pre- 
dicted to occur when 7shO" u £~ 1 ~ 1, with £ cr ~ 0.1 the fraction of 
shock energy carried by suprathermal particles. For 7b — 300, the 
limiting magnetization becomes <r c — 10 -6 , to be compared with 
a u = 5 x 10 -6 B^_ 5 nZ\ in terms of the above fiducial values. 

If o"u 2> cr c , microturbulence cannot be excited at the shock 
front, hence the downstream plasma must be permeated with the 
background shock compressed field only. However, one should not 
expect particle acceleration to proceed and indeed, the PIC simula- 
tions of Sironi & Spitkovsky (2011) indicate no particle accelera- 
tion for a u > 10 -4 , 7b = 15. In order to obtain particle acceler- 
ation, one needs to find a source of turbulence in the downstream 
plasma, that would be able to scatter particles back to the shock 
front faster than they are advected. This possibility cannot be dis- 
carded at the present time but it leads to some form of paradox: 
given that accelerated particles probe a length scale ~ tl behind 
the shock if they interact with a magnetic field coherent over scales 
2> t"l (Lemoine & Revenu 2006, Lemoine et al. 2006), the insta- 
bility needs to produce large scale modes on spatial scales > tl 
with a growth timescale < tl / c. The most reasonable scenario is 
thus to assume that microturbulence exists behind the shock, which 
restricts the possible values of the upstream magnetization. Finally, 
if a u > <j c at some early time, but a u < a c at some late time, be- 
cause <j c oc 7 b ~ 2 while a u remains constant, one should observe a 
very peculiar signature in the light curve due to the initial absence 



of an extended shock accelerated powerlaw (Lemoine & Pelletier 
2011b). 



4 CONCLUSIONS 

Current understanding of the formation of weakly magnetized col- 
lisionless shocks implies the formation of an extended microturbu- 
lent layer, both upstream and downstream of the shock. It is ex- 
pected that the microturbulence behind the shock decays away on 
some hundreds of skin depth scales, and recent high performance 
particle in cell simulations have confirmed this picture. Borrowing 
from these simulations, Section [2] has described the decay of the 
microturbulence as a powerlaw in time, with an index that can take 
mild or pronounced negative values, depending on the breadth of 
the spectrum of magnetic perturbations seeded immediately behind 
the shock. The synchrotron signal of a shock accelerated power- 
law of particles radiating in this evolving microturbulence has been 
calculated, and the results are reported in App. [A] This calcula- 
tion reveals a rather large diversity of possible signatures in the 
spectro-temporal evolution of the synchrotron flux of a decelerat- 
ing relativistic blast wave, commonly parameterized under the form 
F v oc t~ a v~^ . The diversity of signals is associated with the pos- 
sible values of the ratio £dyn/£/i-, which characterizes whether the 
blast has had time to relax to the background shock compressed 
field or not by the comoving dynamical time td yn , as well as the 
different cooling regimes (fast cooling or slow cooling on timescale 
tdyn) and the influence or not of inverse Compton losses on the 
cooling history of particles. One point on which emphasis should be 
put, is that these signatures are potentially visible in early follow-up 
observations of gamma-ray bursts on a wide range of frequencies. 
The detailed synchrotron shapes and indices a, p are discussed in 
App.|X] 

The microturbulence controls the scattering process of par- 
ticles during the acceleration stage and it therefore controls the 
maximal energy of shock accelerated particles. Detailed estimates 
for this maximal energy and for the corresponding energy of syn- 
chrotron photons have been provided. It has been argued that the 
observation of GeV synchrotron photons from gamma-ray burst 
afterglows implies either that the microturbulence decays rather 
slowly (as actually observed so far in PIC simulations), or that the 
spatial extent of the undecayed layer of microturbulence is large 
enough to accomodate both the scattering and the cooling of these 
particles. In the former case, the stretching of the coherence length 
of the microturbulence that accompanies the erosion of magnetic 
power improves the acceleration efficiency and brings it close to a 
Bohm scaling of the scattering frequency. 

Finally, this paper has provided a concrete basis for the phe- 
nomenological models of Barniol-Duran & Kumar (2009, 2010, 
2011a), which interpret the extended GeV emission of a fraction 
of gamma-ray bursts as synchrotron radiation of shock accelerated 
particles in the background shock compressed field. Microturbu- 
lence must play an important role in these models for several rea- 
sons: (1) it is expected to exist behind the shock front of a weakly 
magnetized relativistic shock wave; (2) it ensures the development 
of Fermi acceleration, which cannot proceed without efficient scat- 
tering; (3) it controls the production of photons of energy as high 
as a few GeV. Using the predictions of the spectro-temporal depen- 
dence of F v with a dynamical microturbulence, one can reproduce 
satisfactorily the phenomenology of the models of Barniol-Duran 
& Kumar (2009, 2010, 201 la), as examplified here with the partic- 
ular case of GRB090510. In the present construction, the decaying 
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microturbulent layer permits the generation of high energy photons 
close to the shock front, while the low energy photons are produced 
away from the shock, where the microturbulence has relaxed to (or 
close to) the background shock compressed field. In this setting, the 
broadband early follow-up observations of gamma-ray bursts after- 
glows open an exceptional window on the physics of collisionless 
shocks. More work is certainly needed to study the phenomenology 
of evolving microturbulence and its relation with gamma-ray bursts 
light curves, in particular to define a proper observational strategy 
capable of inferring the characteristics of the microturbulence. 
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The time integral in Eq. lAll integrates the synchrotron power 
along the particle trajectory in the blast, i.e. from shock entry until 
the particle reaches the back of the blast. As it ignores the spatial 
profile of the blast, the upper bound £<jyn is defined up to a factor 
of unity. Such prefactors of order unity are ignored in the following 
calculations, which furthermore always approximate down to bro- 
ken powerlaws. A correct overall normalization is provided at the 
end of the calculation. 

The (comoving) dynamical timescale is defined following 
Panaitescu & Kumar (2000), 
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The spectral power radiated by an electron at time t, of 
Lorentz factor j e is approximated as 
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APPENDIX A: SYNCHROTRON SPECTRA WITH TIME 
DECAYING MICROTURBULENCE 

To compute the synchrotron spectral power of a relativistic blast 
wave, one usually solves a steady state transport equation for the 
blast integrated electron distribution and folds this distribution with 
the individual synchrotron power (Sari et al. 1998). This method 
may be generalized to the present case by including an explicit 
spatial transport term and solving the transport equation through 
the method of characteristics. It appears however more convenient 
to proceed with an approximation of the calculation of Gruzinov 
& Waxman (1999), which integrates over the cooling history of 
each electron in the blast. The present approximation neglects the 
integral over the angular coordinates (Id hypothesis), it neglects 
the hydrodynamical profile of the blast (as in Sari et al. 1998) and 
it neglects the time evolution of the blast physical characteristics 
over the cooling history of a freshly accelerated electron. This re- 
duces the integral to a more tractable form, with a closed expres- 
sion, at the price of reasonable approximations. The angular part of 
the integral indeed leads to a smoothing of the temporal profile of 
emission which may be re-introduced at a later time, see Panaitescu 
& Kumar (2000). The blast evolves on hydrodynamical timescales 
that are longer, or much longer than the radiative timescales. 

The synchrotron spectral power of the blast is written as dis- 
cussed in Section|2T2l 



dN £ 
d7e,0 



^ Ai dE syn 
dt - 



dudt 



(Al) 



where dA r e /d7 e ,o denotes the differential per Lorentz factor inter- 
val of the number of electrons swept by the shock wave by unit 
time. It is defined in Eqs. l 161171 For commodity, the spectral power 
P v is written in the circumbust medium rest frame at redshift z, 
whence the beaming factor 47 b /3, and all frequencies are writ- 
ten in the observer rest frame. This means in particular, that the 
integrands are defined in the comoving blast frame, except dv of 
course, which is an observer frame quantity. The received spectral 
flux F v then reads 



F v = 



(A2) 



with Dl the luminosity distance to the blast. 



Ve = Vp [ 7 eJ SB(t)] 



(A5) 



The peak frequency v p for Lorentz factor 7 e and magnetic field SB 
is defined as (e.g. Wijers & Galama 1999) 



v p [ye; SB] 



3x v e SB 



-7 



47b 



27T m e c 3(1 + z) 



(A6) 



with x p ~ 0.29. For most cases, one might further approximate 
Eq. iA4l as a delta function peaked around v e , but the low energy 
part in i/ 1 / 3 does actually play a role in several specific limits. The 
additional factor of 4/3 in Eq. dA4t ensures proper normalization 
to the synchrotron energy loss rate integrated over frequency. 

The physics of diffusive synchrotron radiation has received a 
lot of attention lately, notably because it might lead to distortions of 
the above spectral shape below and above v v (see e.g. Medvedev 
2000, Fleishman & Urtiev 2010, Kirk & Reville 2010, Reville & 
Kirk 2010, Mao & Wang 201 1, Medvedev et al. 201 1). The magni- 
tude of such distortions is characterized by the wiggler parameter 



eSBXgB 
m e c 2 



(A7) 



If a > 7e, the particle suffers a deflection of order unity in each cell 
of coherence of the turbulent field, hence the standard synchrotron 
regime applies. At the opposite extreme, if a < 1, the deflection 
per coherence cell does not exceed the emission angle of 1 / 7e . This 
leads to diffusive synchrotron radiation, with different slopes at low 
and high frequencies. In the intermediate regime, 1 < a < 7 e , the 
spectrum remains synchrotron like with some departures at v 3> v P 
and v < a^ z v p (Medvedev et al. 2011). 

In the present scenario, one can define a quantity a M immedi- 
ately behind the shock front, in terms of SB^ and A M , 



V8< 



1/2 



A 



— 7b ^ « 20 7m 4 /2 -2^,-o.3Vi , (A8) 



S/Wpi 



with A M i i 

W 4 ElJ/nZl /8 t- 3/8 e 



A M /(10c/^ pi ) and -y m 



3/8 



the minimum Lorentz fac- 



0.3^+, 0.3 

tor of the electron population. For generic parameters, one thus 
finds a M > 7 e or 1 -C a M < y e , but in any case, diffusive 
synchrotron effects can be neglected close to a highly relativistic 
shock front. 

As the turbulence decays, its coherence length evolves, hence 
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a evolves in a non-trivial way, starting from a M . By the time t^- 
i.e. when the turbulence has relaxed to Bd, one finds 



with 



B d 
SB U 



(A9) 



With at < 0, one can check that this value remains much 
larger than unity for generic GRB blast wave characteristics. 
Considering for instance the example at = —3 (fast de- 
cay of the turbulence), a a = 2, one finds a(t^_) ~ 
1.3 x 10 5 B 2 ^ El^t^ 3/S nZf /2A e-^X^^l^. Diffusive syn- 
chrotron effects become important at late times when the blast wave 
has become moderately relativistic, but for the present purposes 
they can be neglected; this justifies the above choice Eq. IA4I for 
dE/dudt. 

The influence of synchrotron self-absorption is discussed 
briefly in Sec. |A5l It does not affect the synchrotron spectra around 
and above the optical range but it may lead to interesting time sig- 
natures at frequencies below a GHz. 



Al 



Gradual decay: 
cooling 



■1 < at < 0; no Inverse Compton 



If —1 < a t < 0, the particle cools gradually in the decaying mi- 
croturbulent field, see Eq. [13] This Section ignores the possibility 
of inverse Compton losses, the effects of which are discussed in 
Sec. |A3l One then defines the critical frequencies 



(A10) 



which is associated to particles of Lorentz factor 7 M + and magnetic 
strength 5B^, and its counterpart 



u p [7 M -;Bd] 



(All) 



The cooling Lorentz factor enters the calculation through the 
upper bound on the integral defined in Eq. dAlt : 7 C denotes as usual 
that which allows cooling on a timescale td yn . To calculate -y c , one 
first define the usual j c ,sb m = 7e t syn [j e ; 8B^\ /t dyn , which cor- 
responds to the cooling Lorentz factor in a homogeneous turbu- 
lence of strength <5B M . The generic notation t syn [y B ; SB] refers to 
the synchrotron loss time for a particle of Lorentz factor 7 C in a 
magnetic field of strength SB. 

Then, if tdyn < t^+, cooling indeed occurs inside the un- 



decayed part of the microturbulence so that 7 C 



7c,5i 



If 



tdyn ^> tn 



cooling instead occurs in the background field Bd, 



meaning 7 C = r y Cy sB l _ l 5B tl / B d . In the intermediate regime t M + < 
tdyn < tjj,-, cooling takes place in the decaying part of the tur- 
bulence. The cooling Lorentz factor is then defined as the solution 
of 



which leads to 



To summarize, 



7c = < 



tsyn [7c ! SB (tdyn)] = tdy 



14-at —at 

7c = 7c,«s fl 7 /J + 



7c, SB^ if tdyn < t M + 

7c,«b' 7/7+ ' if *m+ < *dyn < i M - , 

7c,«S M -gj- if * M - < tdyn 



The characteristic cooling frequency 

v c = Vp [lc\ SB-, C ] , 



(A12) 



(A13) 



(A14) 



(A15) 



SB-, 



max < B d , SB, 



( tdyn\ 



at/2 



(A 16) 



The field 8B lc coiTesponds to the field strength at the point at 
which particles of Lorentz factor 7 C cool through synchrotron radi- 
ation. The latter expression means that if tdyn < tp- , SB lc is given 
by the value of SB at time tdyn, while if tdyn <t fl -, SB Jc — B d - 



The factor (t dyn /t^ + 
with 



|«t/2 



can also be written as (7 C /7 M +) 



St = 



a t 



(A17) 



Note that the value SB-, 



2(1 + at) ' 
can also be understood as the smallest 



value of the magnetic field in the blast. 

The characteristic frequency v m associated to particles of 
Lorentz factor 7 m can also be derived similarly: 



with 



v m = Up [7 m ; SB Jnl ] , 



SB-, — max ^ SB^ C , SB„ { — — 
\7m+ 



(A 18) 



(A19) 



In short, this implies that SB- Ilrl is determined by SB-,,, if j m < 7c 
and by 8B^ (7m/7/ 1 +) ' 5t otherwise, if the 7 m particles can cool 
in the decaying layer. To understand the former value, one should 
note that for at > — 1, most of the power of non-cooling particles 
(7m < 7c) is generated at the back of the blast, since the time 
integrated power oc J idyn dtSB 2 (t) xSB 2 c . 

Depending on how u c , u m , u, L - and z/ M + are ordered one rela- 
tively to each other, one obtains different synchrotron spectra. The 
various possible shapes of these spectra with their characteristic 
indices are summarized in the five cases depicted in Fig. lAll and 
the spectro-temporal indices of F v oc t~ c V~ ti are provided in Ta- 
ble lAll This table assumes a decelerating blast wave in an external 
density profile oc r~ k . 

The following describes in some detail how the spectra are 
obtained, starting from the single particle spectra. Section lA3l dis- 
cusses similar cases under the assumption of strong inverse Comp- 
ton losses. 



Al.l Fast cooling 

Consider first the simplest case of fast cooling, for which u c is 
smaller than all other critical frequencies. The quantity dP v /dN e , 
see Eq. IA1I represents the integral over the history of one cooling 
particle of initial Lorentz factor 7 c ,o- 

If 7e,o > 7m+ , the particle starts to cool in the undecayed 
part of the microturbulent layer where fj,(t) < 1 and continues 
cooling in the decaying part, see Eq.[l3] Straightforward integration 



of Eq.|AT]then leads to u dP v /dN, 



e OC U 



with 



K = < 



'+1/2 if i/ M+ < v < v r ,o,sb„ 

(1 + £>t/2)/( 2 - S t ) if max (i/ M _ ,u c ) <u < u^ + 
+ 1/2 if u c < u < 

1+4/3 if*/<^ c 

(A20) 

The frequency 

Vefi,SB^ = Up [7 e ,o; SB^\ . 
Of course, the spectral power vanishes above i/ e ,o,«B„ ■ 



(A21) 
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Table Al. Spectral (/3) and temporal (a) indices of F v oc t~ a u~@ for various orderings of the 
characteristic frequencies, assuming a decaying microturbulence with — 1 < at < 0, neglect- 
ing inverse Compton losses. The different spectra match those depicted in Fig. lAll Case 1: slow 
cooling, <dy n < tji— \ case 2: slow cooling, t t j Lyn > tp— ; case 3: fast cooling scenario with 
Va— < v c < v m ; case 4: fast cooling scenario with v c < Va— < v m ; case 5: fast cooling 
scenario with u c < u m < v^— . Note that the exact values of the characteristic frequencies vary 
from case to case, see the accompanying text for details. For all cases, —a = (2 — 3p)/4 and 
— f) = — p/2 if u > The quantity k refers to the external density profile n oc r~ k . 



Case 


Frequency range 
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— Q 


Case 1 


v < v m 

U m < V < U c 
U c < U < Vy, + 


1 

3 

1-p 




6(-2+fc) + (-5+fc)c« t 
6(-4+fc) 

fc(10-6p)+24(-l+p) + (-5+fc)(l+p)a t 




2 

-2p+(l-2p)a t 


-2( 


-4+fc)( 


8(-4+fc) 
-2+3p) + (-6+fc(4-5p) + 15p)a t 




4+3a t 






2(-4+fc)(4+3a t ) 


Case 2 


u Uiyi 

Urn < U < Uc 
U c < U < f M _ 

< 1/ < I/ M + 


1 

3 

1-p 
2 

P 
2 

-2p+(l-2p)ct t 


-2( 


-4+fc)( 


2(-3+fc) 
3(-4+fc) 
(-3+fc)(-l+p) 
-4+fc 
-2+fc+3p-fcp 
-4+fc 

-2+3p) + (-6+fc(4-5p) + 15p)c« t 




4+3« t 






2(-4+fc)(4+3c !t ) 


Case 3 


f < f c 

< f < V m 
I'm < f < !>+ 


1 

3 

2+a t 






-4+6fc+3(-5+fc)a t 
6(-4+fc) 
-8+2fc+(-9+fc)a t 




4+3a t 
-2p+(l-2p)a t 


-2( 


-4+fc)( 


2(-4+fc)(4+3a t ) 
-2+3p) + (-6+fc(4-5p) + 15p)a t 




4+3a t 






2(-4+fc)(4+3c- t ) 


Case 4 


f < U c 

U c < U < V M _ 
Up- < 1/ < u m 
V m < U < + 


1 

:S 
1 
2 

2+Q t 






2 

12-3fc 
1 

-4+fc 

-8+2fc+(-9+fc)a t 




4+3a t 
-2p+(l-2p) Qt 


-2( 


-4+fc)( 


2(-4+fc)(4+3a t ) 
-2+3p) + (-6+fc(4-5p) + 15p)a t 




4+3c t 






2(-4+fc)(4+3a t ) 


Case 5 


f < fc 

i^c < U < U m 

I'm < f < !>- 
Vp- < u < u^ + 


1 

:! 
1 
2 
P 
2 

-2p+(l-2p)a t 


-2( 


-4+fc)( 


2 

12-3fc 
1 

-4+fc 
-2+fc+3p-fcp 
-4+fc 

-2+3p) + (-6+fc(4-5p) + 15p)c« t 




4+3a t 






2(-4+fc)(4+3a t ) 



If 7e,o < 7p+, which appears much more likely as discussed 
in Sec. 12. 21 but 7 e ,o > 7 M - (and 7 e ,o > 7c), then u £ fi is obtained 
by solving first for the time t e .o at which the particle actually cools, 
as for Vc and v m : 

tsyn [7e,0i SB (t e ,o)] = *e,0 , (A22) 

which leads to a shifted peak frequency 

Ve.O = V e , . SB ( ) ■ (A23) 

V7m+/ 

Then 

'2 + 2/a t if v e , < v < ^ e ,o,Ss„ 

, (1 + c5 t /2)(2 - 5 t ) ifmax.{vu,-,v c )<v<v e0 

1 — Ou = < 

+ 1/2 if Uc < v < Vp- 

.+4/3 iiv<v c 

(A24) 

For ^ e> o < v < 1^,0, <5B„, the standard index +1/2 has become 
2 + 2/at, which may take large or small negative values depending 
on whether a t lies close to or to — 1 . This index corresponds to the 
radiation of a "non cooling" particle in a changing microturbulent 
magnetic field. 

Finally, if 7 e ,o <S 7 M _ (but 7 e ,o 3> 7 C ), the particle does 
not cool in the decaying microturbulence. The peak frequency for 
v dPv/dNe has moved to 

Ve,0 = fe,o,is„ -^p- , (A25) 



as expected, and for u e fi < v < i/ ei o,5fl„, one finds the slope 
2 + 2/at. Of course, for u c < u < v e ,o, one obtains +1/2. 

Folding the previous results over the particle distribution of 
initial Lorentz factors is straightforward, albeit somewhat tedious. 
For the fast cooling scenario considered here, meaning u c < v m , 
this leads to three generic spectra depicted as cases 3, 4 and 5 
in Fig. IA1I depending on the ordering of v c and i/ m relatively to 
Vfx— '• Vim- < v c < Vm (case 3), u c < Vp— < v m (case 4), 
v c < Vm < Vft- (case 5). This figure does not consider the unlikely 
cases associated to the possibility v m > v^+\ these are briefly ad- 
dressed further below. The distinctive features of these spectra can 
be summarized as follows. 

If Vy,- < Vc, the slope of vF v becomes +1/(2 — St) for 
v c < v < Um', if v c < Vp- however, the slope is +1/2 for u c < 
v < Vp- (fast cooling in B&) and +1/(2 — St) for v^- < v < v m 
(fast cooling in decaying turbulence). For v m < v < v^+, the slope 
is (2 — p)/(2 — St). The general trend of decaying microturbulence 
is to produce flatter synchrotron spectra than in a homogeneous 
magnetic field, due to the stretch in frequency associated to cooling 
in regions of different magnetic field strengths. 

The maximal frequency c maj( associated to 7 max is calcu- 
lated in a similar fashion to v m . This frequency does not appear 
in Fig. lAll for the sake of clarity but its impact can be described 
as follows. If f m „ > i7i+, spectral power vanishes above i/ mlx . 
Otherwise, uF v has spectral index 2 + 2/at for f m „ < v < 
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Figure Al. Generic synchrotron spectra in time decaying microturbulence 
with — 1 < at < 0, neglecting inverse Compton losses, with the spec- 
tral indices as indicated. Case 1: slow cooling scenario (7 C > 7 m ) with 
tdyn < tpt— > m which case the turbulence has not had time to relax to the 
background value Bd', case 2: slow cooling scenario with t^ vn > £„_ ; 
case 3: fast cooling scenario with Va— < v c < i/ m ; case 4: fast cool- 
ing scenario with u c < < !^ m ; case 5: fast cooling scenario with 
v c < I'm < i^u— ■ See the discussion in Sec. lAll for the precise definitions 
of the various frequencies and for modifications to these spectra if u c or v m 
exceeds or if i/ maJ[ is taken into account. 



Vp [7max; t>B M ], which as before may take values close to zero if 
at is close to —1, and it vanishes beyond v p [7ma X ; SB^]. 
The peak power vF v is radiated at v m as usual, with 

0.28 4 



vF v \ 



4%D 2 L 3 



(A26) 



The above relates the maximum power to the incoming electron 
energy per unit time, which avoids specifying the value of the mag- 
netic field in which cooling takes place. The numerical prefactor 
0.28 matches the prefactors derived in Panaitescu & Kumar (2000) 
for non-decaying turbulence, Lorentz beaming is included through 
the factor 4/y^/3 and N e has been defined in Eq. d!6t . 

A1.2 Slow cooling 

For high energy particles with 7 ej o > 7c the individual spectra of 
the quantity v dP„ /dN e mimic those discussed in the fast cooling 
section before and this discussion thus concentrates on the bulk of 
electrons for which 7^0 < 7c- Such electrons do not cool substan- 
tially anywhere in this slow cooling limit, therefore one finds either 



a slope +4/3 at low frequencies, or a slope 2 + 2 /at at high fre- 
quencies corresponding to the changing magnetic field. There are 
two frequencies associated to ^/e,a'- Vefi.SB^ as before, correspond- 
ing to the peak frequency of emission when the particle experi- 
ences the undecayed 5B M , and f e ,o the frequency to be calculated 
in the lowest magnetic field found downstream, i.e. at the back of 
the blast. Then 



,0 = v v [7e,o; SB*y 



(A27) 



provided 7 e ,o < 7 M +. 

As before, folding over the particle initial Lorentz distribution 
leads to the all-particle spectra. The possible spectra are displayed 
as cases 1 and 2 in Fig. lAll depending on the ordering of td yn vs 

case 1 if t<x yn < case 2 if tdyn > t/j-, which implies 
v c < Vfi- and 8B lc — Bd- If v c > v^+, one would of course 
recover the spectrum of a standard slow cooling scenario in a ho- 
mogeneous turbulence of strength 5B M . 

The distinctive features of the vF v spectra can be summa- 
rized as follows: the slope 1 — p/2 of the fast cooling part of 
the particle population has been turned into (2 — f>)/(2 — 5t) for 
max (y c , < v < V/J.+. The +4/3 slope for v < ^ m remains 
unchanged, just as (3 — p)/2 for u m < v < v c , or 1 — p/2 for 

< v (assuming < ^ max of course). 

The peak power is now radiated at v c , with 



vFJ 



0.28 4 2 ■ 2 

7bA'e7cW- e C 



47rD| 3 
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A2 Rapid decay: a t < 1; no inverse Compton cooling 

This Section now considers the limit in which the energy density 
stored in the microturbulence decreases faster than t . Inverse 
Compton losses are neglected here; their impact is discussed in 
Sec. IA4I The crucial difference between the limit a t < —1 and 
that discussed in the previous Section has to do with the cooling 
history of a particle. In the present case, either the initial Lorentz 
factor y e fi > 7 M +, in which case the particle cools down to 7 M + 
after crossing the undecayed part of the microturbulent layer, or 
7e,o ^ 7m+> in which case it does not cool anywhere in the decay- 
ing microturbulent layer, see Eq. (1131 . In this latter case, the par- 
ticle eventually cools in the background magnetic field, provided 

^dyn ^ "11— • 

The cooling Lorentz factor and its corresponding frequency 
should therefore be defined as follows. If tdyn < one recov- 
ers a trivial case as it means that the turbulence has not relaxed 
beyond <5£? M , hence the turbulence is homogeneous downstream. 
This case is not discussed further here. If t M _ < td yn , the tur- 
bulence has relaxed to B^ by td yn ; in this case, j c is defined as 
usual in terms of tdyn and Bd- Indeed, cooling cannot take place 
in the decaying part of the microturbulent layer, but cooling is pos- 
sible in the background shock compressed field Bd- The cooling 
frequency then reads v c — v p [j c ; Bd]. In the intermediate limit, 
tft+ < tdyn < t^-, the turbulence has not had time to relax down 
to Bd. The cooling Lorentz factor remains undefined; however, 
one can understand the spectra obtained with the correspondence 
7c — > 7/j+ and v c — > since particles above 7 M+ do cool down 
to 7, 1+ in 5B M . 

For a particle of initial Lorentz factor 7 e ,o, one should define 
two critical frequencies: f e ,o,5s„ = v P [7e,o; <5B M ] and v e ,o.B d = 
v p [7e,o; Bd]. Note that the Lorentz factors 7 M + and 7 M _ and their 
corresponding frequencies and f M _ remain unchanged here. It 
proves necessary to define a new frequency associated to particles 
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of Lorentz factor radiating in the lowest magnetic field SB-, 
with as before 8B-, 



One also define 



Bd if tdyn > tp-'. 

v^o = u p [7, 1+ ; SB lc ] 

Vm,6B„ = fp [7m! SB U ] 



(A29) 



(A30) 



The expression of the characteristic frequency i/ m is given further 
below, case by case. 



Spectra v AP V / 'dN £ oc 



integrated over the cooling his- 



tory of a particle of initial Lorentz factor 7 e ,o show: 



l-b v 



f +1/2 
+ 1/2 

min(2 + 2/a t , 4/3) 



if V,t+ < V < f e ,o,SB p 

if v c < v < min (u u -, v e ,o,B d ) 
if min (v u -, v e ,o,B d ) < v 
and v < min (vy,+, v^q^b^) 
(A31) 

To understand the latter slope, one may recall that in a decaying mi- 
croturbulence, radiation in a region of small extent but high mag- 
netic power competes with radiation in a region of large extent at 
small magnetic power. If a t < —3, decay is so fast that most of the 
radiation is produced in 8B U and one collects at low frequencies 
the +4/3 tail. If however, —3 < at < — 1, the radiation produced 
by the particle as it crosses the decaying part of the turbulence dom- 
inates this tail and the slope becomes 2 + 2/a t . This latter can be 
much flatter, possibly giving rise to a flat energy spectrum in the 
limit at — > — 1. 

After folding over the particle population, one obtains the 
generic spectra depicted in Fig. IA2l for —3 < a t < — 4/(p + 1) 
and in Fig. lA3l for — 4/(p + 1) < at < — 1. The salient features of 
these spectra can be summarized as follows. 

Consider for simplicity the limit —3 < at < — 4/(p + 1). 
The other limit — 4/(p + 1) < at < —1 can be understood in a 
similar way, while the limit at < —3 follows from the former after 
replacing 2 + 2/at with 4/3. If fdyn < t u -, the turbulence has not 
had time to relax down to Bd- At high frequencies > ^ m ,iB , the 
spectrum then takes the form of a slow cooling scenario in a homo- 
geneous turbulence of strength SB^, with v c — > As discussed 
above, the turbulence decays so fast that the early emission in the 
region of high magnetic power dominates that further away from 
the shock. At frequencies t> m < v < !> m ,,5B , one collects the low 
energy extension with slope 2 + 2/at, instead of 4/3 as discussed 
before. Note that the frequency ^ m = u p [y m ; 5B lc ]. For v < v m , 
one recovers the slope +4/3 as expected. This case is denoted case 
1 in Fig.lA2l 

In case 2, one now assumes fdyn > t a - (relaxed turbulence) 
and u c > v m (slow cooling). The afterglow comprises two con- 
tributions: one associated to slow cooling in SB^, as before; plus 
a second one associated to slow cooling in the background Bd- 
This latter is indicated in dashed lines in Fig. IA2I it exhibits a 
cut-off above v u o, since there are no particles with Lorentz fac- 
tor above j u + beyond t u +. If 7 max < 7 M +, the cut-off would 
of course take place at v p [7maxj Bd] - The characteristic frequency 
ZAn = v p [7 m ; Bd]- 

The peak power for the vF v .b a component associated to cool- 
ing in Bd is standard, while the power uF Vi gBu, related to the de- 
caying microturbulence can be written at f M -|_ as 



0-28 4 2 .-. 

- Jb N el , + m e c 



.2 ( 7m+ 

V 7m 





Figure A2. Generic synchrotron spectra in time decaying microturbulence 
with —3 < at < — 4/(p + 1), neglecting inverse Compton losses, with 
the spectral indices as indicated. Case 1: scenario with id yn < t^—, in 
which case the turbulence has not had time to relax to the background value 
Bd', case 2: slow cooling scenario with td yn > ty,— \ case 3: fast cooling 



scenario with td y 



> t u 



See the discussion in Sec. IA2l for the precise 



definitions of the various frequencies. To obtain the spectra for at < —3, 
it suffices to carry out the replacement 2 + 2/at — > 4/3. The dashed line 
represents the secondary synchrotron component associated to cooling in 
the background shock compressed field Bd, whenever tdyn > tfi— ■ 



The ratio of energy fluxes at their respective peaks reads 

vF u ,b a ( 7c = {SB't,, 

7aH 



(A33) 



vF v .SB^ V7m+7 V ^dyn 

and it scales as one would expect with the ratio of the product of 
synchrotron power times the exposure to the magnetic field. 
One may also calculate the energy flux ratio at i/ m : 

ZAn-f^m ,B d fdyn 

l / m-ti/„.i5B.. t u - 



(A34) 



In case 3, one now assumes a similar configuration with 
fdyn > t u - (relaxed turbulence), but u c < v m , meaning fast cool- 
ing in Bd- The only difference with the previous spectrum corre- 
sponds to that change of cooling regime for the component vF v ^b a ■ 
As discussed before, the bulk of electrons cannot cool in the micro- 
turbulent field (as long as 7 m < 7 M +). 

Figure IA3I shows the corresponding spectra for the case 
— 4/(p + 1) < at < — 1. The discussion is very similar to the 
previous one for at < — 4/(p + 1) and the differences are as fol- 
lows. For — 4/(p + 1) < at < — 1, the spectral index 2 + 2/at 
is softer than (3 — p)/2; the low energy tail of index 2 + 2/at of 
the high energy population thus extends from down to v u o as 
it dominates the slow cooling contribution of the bulk of electrons. 



(A32) 



The spectro-temporal indices of F v oc f 
values of a t are provided in Table IA21 



for these two 



18 M. Lemoine 



Table A2. Spectral (/3) and temporal (a) indices of F v oc t~ a v~P assuming a decaying microturbulence 
with —3 < at < — 1, with negligible inverse Compton losses. In cases 2 and 3, one must superimpose 
a synchrotron component associated to cooling in the background shock compressed magnetic field; the 
spectro-temporal slopes given here concern only the synchrotron component associated to the decaying 
microturbulent layer, not the latter. Case 1: id yn < t/j— ! case 2: id yn > t/j— with slow cooling in 
the background shock compressed field; case 3: t^ yn > t M _ with fast cooling in the background shock 
compressed field. The corresponding synchrotron spectra are shown in Fig. lA2l for — 3 < at < 4/(p + 1) 
and in Fig. IA3l for — 4/(p + 1) < at < — 1. For all cases, —a = (2 — 3p)/4 and — /3 = — p/2 if 
" > The quantity k refers to the external density profile n oc r~ k . 

Case Frequency range —/J — a 



Case 1 ^ < f n 



1 6(-2+fc) + (-5+fc)c t 

3 6(-4+fc) 

-4/(p+l)<at<-l] . m <.<^o V Mio-6p) + 24(-i+pH.(- 5+fc )(i +rtat 

,. „ ^ ,. ^ ,. , 1,2 2(6+fc) + (2-fc(-2+p)+9p)at 



Cases 2 and 3 v < f m 



1 -6fc + (3 + fc)at 

3 6(-4+fc)a t 

[-4/(p+l)<a t <-l] ^ m <^<^o ±? ' 2H( 2(-4 + »; aP '°' 

,y „ <f u u , i i 2 2(6+fc) + (2-fc(-2+p)+9p)a t 

Case 1 „ <" „ 1 6(-2+fc) + (-5+fc)a t 

Lase i K < l/ m 3 6(-4+fc) 

r o ^ a//- i 1\1 ^ ^ 1 i 2 6(-4+fc) + (-7+3fc)a t 
[-3 < a t < -4/(p + 1)] u m < v < v m> s Bll 1 + - 



"m,iB„ < " < "m + 



at 2(-4+fc)a t 
1-p 2-3fc-12p+3fcp 
2 16-4fc 



1 -6fc + (3 + fc)a t 

3 6(-4+fc)a t 

r a ^ ^ a If i 1Y1 i i 2 6(-4+fc) + (-7+3fc)a t 

[-3 < at < -4/(p + 1)] i> m < i/ < f m , 5S(J 1 + — 2{U+k) at 



Cases 2 and 3 v < v 

I'm < 

u ui,SB u < V < + 



Qt 2(-4+fc)c* t 
1-p 2-3fc-12p+3fcp 
2 16-4fc 





Figure A3. Same as Fig.lA2lfor -4/(p + 1) < o t < -1. 



A3 Gradual decay at > — 4/(p + 1), with dominant inverse 
Compton cooling 

This Section assumes —4/ (p + 1) < at < and it assumes that 
the Compton parameter Y 3> 1 for all Lorentz factors, everywhere 
in the blast. One defines 



Y a = 



Ur, 



" " <5B2/(8Tr) ' 
The cooling time of a particle is then defined as 

1 



tcool(7e,o) 



1 + Y, 



■tsyn [je,0\ SBfj 



(A35) 



(A36) 



and it is homogeneous throughout the blast. This simplifies the 
cooling history of a particle: 



7e,0 if t < t coo \ (7e,o) 

icool (7e,o) 
, 7e '° i 



if t > tcool (7e,0) 



(A37) 



The definitions of the critical Lorentz factors must be adapted 
to this case. One defines a Lorentz factor 



7m+ 



1M 



(A38) 



such that inverse Compton cooling takes place at the end of the 
undecayed microturbulent layer. Similarly, one defines ■y )l - = 
7 M +t M +/t M _. The associated frequency !/ M+ (resp. is defined 
as before in terms of 7 M + (resp. 7 M _). 

The cooling Lorentz factor is defined as 



7c = 7m+ 



£p+ 

^dyn 



(A39) 



Eqs. |A15|A16l for the cooling frequency u c and the definition of 
SB lc remain valid. Regarding v m , Eg. I A 181 remains valid but the 
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definition of SB lm must be modified to account for inverse Comp- 
ton losses: 



SB,, 



max i SB 7 , SB^ ( ) 



-at/2 



(A40) 



The index of v AP v /AN e differs from the standard case of ho- 
mogeneous turbulence only if the particle radiates in the changing 
magnetic field while it is cooling in the radiation field. This applies 
to the spectral domain max [y^-, v c ) < v < min {v e ,o, ^n+), in 
which the index of v dP v /dN e becomes ( 1 — at)/ (2 — at/2). The 
frequency i/ e> o is defined in terms of 7e : o as Vm in terms of 7 m . 

After folding over the particle Lorentz distribution, one ob- 
tains the generic full all-particle spectra represented in Fig, IA4l As 
before, this figure ignores a possible maximal frequency ^ max and 
assumes that v c and u m are smaller than The discussion is 
very similar to that given in Sec. IA1 l and is not be repeated here. The 
significant differences lie in the spectral indices: (2 — p)/(2 — 8t) 
has become (2 — p — at)/ (2 — at/2), while 1/(2 — St) (otherwise 
1/2 for standard fast cooling) has become (1 — at)/ (2 — at/2). 
Inverse Compton losses take away a large fraction of the dissipated 
energy, so that the synchrotron peak power is modified as follows, 
for slow cooling (cases 1 and 2): 



vFJ 



0.28 4 2 ■ 

4W3 7biVeW ( ~ 



i-p 



SB?, 



1 + SB* 

and for fast cooling (cases 3, 4 and 5): 



(A41) 



vFJ 



0-28 4 a . 2 

4^|3 7biVe7mmeC X 

1 SB* 



1 + Yu SBl 



(A42) 



The ratios of magnetic energy densities SB^/SB* (resp. 
SB^^/SB*) that appear in the both expressions, yield the proper 
Y Compton parameter at the location at which most of the cooling 
of particles of Lorentz factor 7 C (resp. y m ) occurs. 

The corresponding spectro-temporal indices for F v oc 
t~ a u~ p are given in Table lA3l 



A4 Rapid decay —3 < a t < — 4/(p + 1), with dominant 
inverse Compton cooling 

Finally, one must consider the possibility that —3 < at < — 4/(p+ 
1) with dominant inverse Compton losses, as in the previous Sec- 
tion |A3l The generic spectra obtained after folding over the particle 
population are depicted in Fig. lA5l The critical frequencies are de- 
fined as in the previous Section |A3l 

The characteristic features of these spectra mix those of 
Sec. [A3] [gradual decay, at > — 4/(p + 1) with inverse Compton 
losses] with those of Sec. IA2I (rapid decay, at < —1, no inverse 
Compton losses). They can be summarized and understood as fol- 
lows. 

Cases 1 and 2 depict a slow cooling scenario 7 m < 7 C with 
either tdyn < t^- (case 1) or tdyn > t^- (case 2). Because the 
bulk of electrons do not actually cool, the spectra obtained match 
those of cases 1 and 2 in Fig. lA2l for which inverse Compton losses 
are absent. 

Cases 3 and 4 depict fast cooling regimes with different order- 
ings of v Q relatively to f M _, but v m > v^-. The peak power still 




Figure A4. Generic synchrotron spectra in time decaying microturbulence 
with — 4/(p + 1) < at < 0, assuming that inverse Compton losses dom- 
inate over synchrotron losses everywhere in the blast. Case 1 : slow cooling 
scenario (_j c > 7 m ) with tdyn < tp— > in which case the turbulence has not 
had time to relax to the background value \ case 2: slow cooling scenario 
with tdyn > t/i— ; case 3: fast cooling scenario with v^— < u c < u m ; case 
4: fast cooling scenario with v c < f„_ < u m ; case 5: fast cooling sce- 
nario with v c < u m i n <u^—. See the discussion in Sec. lA3l for the precise 
definitions of the various frequencies and for modifications to these spectra 
if v c or v m exceeds or if i/ maI is taken into account. Note that the 
spectral slope (2 — p — Qt)/(2 — at/2) takes positive values whenever 
a t < 2 — p. 



occurs at due to the rapid decay of the microturbulence and 
the spectra keep a shape similar to those of cases 1 and 2, up to the 
fast decaying index at low frequencies v c < v < v m . Case 4 dif- 
fers somewhat from its counterpart when inverse Compton losses 
are neglected, which is represented as case 3 in Fig. lA2l because it 
does not reveal an additional synchrotron fast cooling component 
in the background shock compressed magnetic field. This absence 
is directly related to the presence of inverse Compton losses, which 
cut the electron distribution down to 7^_ < j m by the end of the 
microturbulent layer, so that the electron distribution is mostly mo- 
noenergetic at that point and cooling in the background magnetic 
field ensues with generic slope 1 /2. 

In contrast, 7 m < 7 M - in case 5. Then the electron distribution 
maintains a powerlaw shape between 7 m and 7 M _ at the end of the 
microturbulence layer, so that an additional component associated 
to cooling in the background field emerges, as shown in Fig. |A5| 

The corresponding spectro-temporal indices for F v oc t~ a v^ 
are given in TablelA4l 
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Table A3. Spectral (/3) and temporal (a) indices of F v oc t~ a u~P for various orderings 
of the characteristic frequencies, assuming a decaying microturbulence with — 4/(p + 
1) < at, with dominant inverse Compton losses everywhere in the blast. Case 1: slow 
cooling with t^ yn < case 2: slow cooling with td yn > i/j— ! case 3: fast cooling 
scenario with < u c < i/ m ; case 4: fast cooling scenario with u c < Up- < !/ m ; case 
5: fast cooling scenario with u c < i/ m in < f/j— . The corresponding synchrotron spectra 
are shown in Fig. IA4I For all cases, —a = (2 — 3p)/4 and —f$ = —p/2 if u > 
The quantity A; refers to the external density profile nar"'. 



Case Frequency range —/J —a 



Case 1 ;/ < !/ m 


1 

3 


6(-2+fc) + (-5+fc)a t 
6(-4+fc) 


l/m < V < V c 


1-p 

2 


fc(10-6p)+24(-l+p) + (-5+fc)(l+p)a t 
8(-4+fc) 


u c < u < Up + 


2p+a t 


2(-4+fc)(-2+3p) + (2-fc(-2+p)+9p)a t 


-4+a t 


2(-4+fc)(-4+a t ) 



Case 2 i/ < f m 

i/ m < ^ < 

Uc < U < Vfj,- 

V u - < U < v u 



1 2(-3+fc) 

3 3(-4+fc) 

1-P (-3+fc)(-l+p) 

2 -4+fc 

_p 2+(-3+fc)p 

2p+q t 2 2(-4+fc)(-2+3p) + (2-fc(-2+p)"+9p fc )at 

-4+a t 2(-4+fc)(-4+a t ) 



Case 3 v < u c 

u c < u < I/, 



1 -4+6fc+(-5+fc)a t 

3 6(-4+fc) 

2+q t 2(-4+fc) + (ll+fc)n t 

-4+eti 2(-4+fc)(-4+a t ) 

, „ , ,, 2p+ Qt 2(-4+fc)(-2+3p) + (2-fc(-2+p)+9p)a f 

^ " + -4+a t 2(-4+fc)(-4+« t ) 



Case 4 f < u c 

U c < U < Up- 
Up.- < U < Um 
U m < V < Up + 



1 2(-l+fc) 
3 3(-4+fc) 
1 l-h 



2 -4+fc 

2+at 2(-4+fc) + (ll+fc)n t 

-4+Q t 2(-4+fc)(-4+a t ) 

2p+c t 2(-4+fc)(-2+3p) + (2-fc(-2+p)+9p)a t 

-4+ Qt 2 (-4+fc)(-4+a t ) 



Case 5 u < u c 


1 

3 


2(-l+fc) 
3(-4+fc) 


u c < u < l/m 


1 

2 


1 — fc 
-4+fc 


I^m < V < Up- 


p 

2 


2+(-3+fc)p 
-4+fc 


Up- < U < Up + 


2p+c*t 


2(-4+fc)(-2+3p) + (2-fc(-2+p)+9p)a t 


-4+aj 


2(-4+fc)(-4+a t ) 



The peak power for the synchrotron component associated to 
the microturbulent layer occurs at Vp+ in all cases with 

P I 0-28 4 2 - 7 M+ 2 /7m+ N 

(A43) 

The ratio between the peak powers of each synchrotron component 
for case 5 reads 



(A44) 



which resembles Eq. lA44l up to an additional factor B^/5B^ at the 
benefit of the high energy microturbulent component, due to the 
competition between synchrotron and inverse Compton losses. 

A5 Synchrotron self-absorption 

At very low frequencies, the synchrotron spectrum may be modi- 
fied by opacity effects. The absorption break frequency generally 
lies well below the optical domain, e.g. for a homogeneous turbu- 
lence, k = (constant density profile) and a slow cooling regime, 
one finds i/ a b s ~ 3 x 10 7 Hz El^ 5 n 3 /^e 1 J 5 _ 2 e~]_ 3 z^ 1 3 (see 
Panaitescu & Kumar 2000). In this standard case, the frequency in- 
dex of F v is 2 below the absorption frequency, e.g. Granot et al, 
(1999) and references therein. 

Accounting for the decaying microturbulence modifies the sit- 
uation as follows. All calculations are performed in the comoving 



downstream frame, as indicated by the primes. Moreover, the self- 
similar profile of the blast and the secular evolution of the blast 
characteristics are neglected, as in the rest of App. [A] One then de- 
fines the synchrotron self-absorption coefficient (Rybicki & Light- 
man, 1979): 



8ivm 



/e7e 



d ( 1 dn' 



d-Je \7l d 7e 



(A45) 



with P' v i i e = dE/dv'dt' the spectral synchrotron power emit- 
ted per electron per frequency interval, as defined in Eq. IA4l up 
to the change of frame for the frequency. Here however, one con- 
siders only the low frequency part oc v 1 ^ since one is interested 
in frequencies well below i/ m . The electron distribution function 
dn' /d~/ e depends on the distance to the shock front when cooling 
is efficient. Since only the minimum Lorentz factor is of impor- 
tance in the calculations of synchrotron self-absorption, one can 
approximate the distribution as a unique powerlaw, 



dn' 1 ( le y q . 
— = — — n(q-l) 
d7e 7i \Hj 



(A46) 



with 7^ = 7m at distances such that cooling is inefficient for parti- 
cles of Lorentz factor 7 m , and ji position dependent at larger dis- 
tances from the shock front. Of course, in a slow cooling regime, 
7< = 7m everywhere in the blast. The index q departs from s only 
when cooling becomes efficient; its exact value does not affect the 
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Table A4. Spectral (/3) and temporal (a) indices of F v oc t~ a v~& for various or- 
derings of the characteristic frequencies, assuming a decaying microturbulence with 
—3 < at < — 4/(p + 1), with dominant inverse Compton losses everywhere in the 
blast. In case 2 and 5, one must superimpose a synchrotron component associated to 
cooling in the background shock compressed magnetic field; the spectro-temporal slopes 
given here concern only the synchrotron component associated to the decaying microtur- 
bulent layer, not the latter. Case 1: slow cooling with t c \ yn < f, M _; case 2: slow cool- 
ing scenario with id yn > tp~\ case 3: fast cooling scenario with u^- < u c < u m ; 
case 4: fast cooling scenario with u c < v^— < Um\ case 5: fast cooling scenario with 
< ^min < ^n— ■ The corresponding synchrotron spectra are shown in Fig. IA5I For 
all cases, —a = (2 — 3p)/4 and — /3 = —p/2 if u > 17,+ . The quantity k refers to the 
external density profile n <x r~ k . 

Case Frequency range — —a 

Case 1 u < u m 

Case 2 u < i/ m 

Vm,SB,_, <V< U^ + 



-4+6fc+(-5 + fc)q t 
6(-4+fc) 
2(-4+fc) + (H+fc)"t 
2(-4+fc)(-4+a t ) 
6(-4+fc) + (-7+3fc)a t 
2(-4+k)a 4 
2-3fc-12p+3fcp 
16-4fc 



Case 4 < 

I/ c < U < V/A— 
Up.- < U < Um 

v m <v < v m j Bli 

V m ,SBp <V < Up + 

Case 5 u < v c 

U c < U < Um 

Um < V < l/ m ,i B(1 

"miflj, < V < Up + 



1 6(-2 + fc) + (-5 + fc)c t 

3 6(-4+k) 

_2_ 6(-4+fc) + (-7+3fc)a t 

a t 2(-4+fc)a t 

1-p 2-3fc-12p+3fcp 

2 16-4*; 



1 -6fc+(3+fc)a t 

3 6(-4+fe)a t 

_2_ 6(-4+fc) + (-7+3fc)a t 

a t 2(-4+fe)a t 

1-p 2-3fc-12p + 3fcp 

2 16-4fc 



Case 3 u < u c 

U c < U < Um 
Um < U < U„ 



1 + 



i+a t 
2 



"m,«B„ < " < "p + 



_L -16 + 6fc-5(-5 + fc)q f 

3 24(-4+fc) 

1 1-fc 



1 + 



2 -A+k 

2 + at 2(-4+fc) + (ll + fc)q t 

-i+a t 2(-4+fc)(-4+«t) 

_2_ 6(-4+fc) + (-7+3fc)a t 

a t 2(-4+fc)a t 

1-p 2-3fc-12p + 3fcp 

2 16-4fe 



1 + 

I 

2 



1 -6fc+(23+fc)a t 

3 6(-4+fc)a t 

1 -2fc+(ll-3fc)at 

"2 2(-4+fc)at 

2_ 6(-4+fc) + (-7+3fc)a t 

at 2(-4+fc)at 

1-p 2-3fc-12p+3fcp 



following up to factors of order unity. The above integral leads to 

with vl = u' p [<5-B;7;]. This local absorption coefficient depends 
on position both through 7; (when cooling is efficient) and through 
SB, in the presence of decaying turbulence. One may thus rewrite 
the absorption coefficient as 

/ c R \ 2/3 / \ -5/3 

(£) ■ (A48 » 

with a'p v , the absorption coefficient defined according to Eq, IA47l 
with SB — > 5 Bp and 7, — > 7 m . This formulation allows to write 
q^/ as a broken power law function of the distance to the shock 
front. 

The emission coefficient j' v , is defined as 

>* = hf (A49) 

assuming isotropic emission, and it can be recast similarly to a' v , 



into 

with ~ 0.4fi 8 *B M (m.c?)- 1 n'(i//i^, m ) 1/8 and 1/^ = 
[5-B M ;7 m ]. 

The specific intensity along a ray path obeys the radiative 
transfer equation 

c\T' , 

*±-± = -cl v ,l' v , + j' v , , (A51) 
as 

with the usual formal solution 

I' v ,{x = 0) = P'' 15 e^' ■ (A52) 

Here dr, y / = a' v ,dx, x > denotes the distance from the shock 
front, and t v i (x = 0) = by convention. The parameter T„/ b con- 
sequently represents the total optical depth of the blast which can 
be derived by integrating the broken power law form of Eq. IA48I 
over x, from x = up to x^ ~ f/7b ~ id yn /c. The above formu- 
lation of the solution of the equation of radiative transfer is particu- 
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Figure A5. Generic synchrotron spectra in time decaying microturbulence 
with —3 < at < — 4/(p + 1), assuming that inverse Compton losses 
dominate throughout the blast. The spectral indices as indicated. Case 1: 
slow cooling scenario with t^ yn < t M _, in which case the turbulence has 
not had time to relax to the background value B d ; case 2: slow cooling 
scenario with td yI1 > i^— ; case 3: fast cooling scenario with v^— < 
Vc < case 4: fast cooling scenario with u c < v^- < Vm \ case 5: 
fast cooling scenario with v c < i/ m < Vu,—. See the discussion in Sec. lA4l 
for the precise definitions of the various frequencies and for modifications to 
these spectra if u c or f m exceeds , or if v ma x is taken into account. The 
dashed line represents the secondary synchrotron component associated to 
cooling in the background shock compressed field as in Fig. lA2l 

larly appealing because it leads to a simple evaluation of I' v , (0) in 
various cases. 

In particular, if the regime is slow cooling, 74 = j m every- 
where and I oi v i becomes uniform in the blast, as in the homo- 
geneous (slow cooling) case. Then one finds 

I',A0) = ^ (l-e- T '^) . (A53) 

In the optically thick limit t v i . b 2> 1, the ratio j' v , /a' v , leads 
to the standard slope 2 for F v below the absotption frequency 
^abs t>> which is defined by the condition 7V jb = 1 at ^ ba b . 
Note that r„/ b shares with the absotption coefficient a' v , the scal- 
ing oc !/'~ 5,/3 . In the optically thin limit T„/ jb <C 1, meaning 
v' 3> I'absjb- one recovers the slope 1/3 since t v i >b oc a' v ,. The 
non trivial dependence of v'^ hsh on time leads to a non trivial time 
dependence of the flux below the break frequency. This dependence 
is not discussed here, but can be calculated through 7V ib . 



The fast cooling regime leads to a more complicated shape of 
the spectrum in the transition region between the asymptotic thin 
and thick regimes (for a similar discussion with a homogeneous 
turbulence, see Granot et al. 2000). One needs to define here an 
intermediate optical depth, 

T v >, m = / ov dx , (A54) 

Jo 

with a- coo i. m the location at which particles of Lorentz factor 7 m 
start to cool. The fast cooling regime implies that by the back of 
the blast, such particles have cooled down to <y c . One now defines 
an absorption frequency v' &hs m such that r v i >m = 1 at v' aha m . 
Because the optical depth is an increasing function of distance 
to the shock front, r v i ?m < t v i ib (possibly r, y /, m <C t v i , b ) and 
^abs.m < ^abs,b- Therefore, at frequencies u' < ^ bs>m , the shell 
is optically thick for x ^ x C oo\,m and the standard synchrotron 
self-absorption spectrum emerges. At frequencies v' > v' aba b , the 
whole shell is optically thin to synchrotron radiation, so that the 
standard oc i/ 1//3 spectrum calculated above applies. However, for 
^abs m < v ' < v 'a!as b> opacity competes with local emission at 
distances > x coo i,m, and the overall spectrum departs from either 
asymptote. The exact spectral shape can be derived from Eq. IA53I 
using the calculated value of T„/ b . 
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